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ABSTRACT 


Reported  here  are  the  results  of  a  three  year  study  of  the  high 
rate  volumetric  Properties  of  Snow.  In  addition  to  the  report  presented 
here,  the  project  resulted  witli  over  ten  publications  (Brown  1179-1.981  , 
Bowles  1981,  Rose,  1181)  which  are  included  in  the  list  of  reference's. 

Tin-  constitutive  behavior  of  dry  snow  is  analyzed  and  described  mathe¬ 
matically  by  means  of  two  microdynamica 1  models.  These  are  then  evaluated 
by  means  of  previous  experimental  data.  A  number  of  applications  are 
then  investigated.  These  applications  include  vehicle  mobility  in 
shallow  and  deep  snowpack,  steady  shockwaves,  and  nonsteady  shockwaves. 
Finally,  an  electromagnetic  stress  wave  generator  was  designed,  constructed 
and  used  in  a  shockwave  testing  program  on  snow.  experimental  results 
arc  presented  and  discussed. 
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INTRODUCTION 


I_. 

This  report  is  concerned  primarily  with  the  high  strain  rate  volumetric 
properties  of  snow.  This  work  was  originally  oriented  toward  the  analysis 
of  shockwaves  in  snow,  but  as  the  project  progressed,  it  became  apparent 
that  the  study  would  have  to  include  a  major  effort  to  evaluate  the 
volumetric  properties  of  snow.  This  property  of  snow  must  be  known 
before  a  number  of  probems  involving  snow  can  be  studied. 

The  mechanical  properties  of  snow  have  been  under  active  investigation 
for  over  forty  years.  Mellor  (1974,  1977)  lias  made  thorough  reviews  of 
past  investigative  work  on  the  properties  of  snow.  Salm  (1981)  has 
recently  made  a  more  up  to  date  review  and  has  described  the  current 
state  of  the  art. 

Marly  work  (1940-1960)  was  directed  at  evaluating  the  linear  properties 

of  snow.  Work  typical  of  this  period  includes  work  by  Kinosita  (1967) . 

Later  work  by  Lee  (1961,1963),  and  Salm  (1974)  extended  the  earlier  work 

on  the  linear  properties  of  snow.  Most  of  this  work  concentrated  on  the 

linear  viscoelastic  properties.  At  low  enough  strain  rates,  say  at 
-7  -1 

rates  less  than  10  s  ,  snow  does  exhibit  linear  properties,  although  the 
deformation  is  generally  inelastic.  Most  linear  constitutive  relations 
used  to  describe  the  behavior  of  snow,  have  the  forms 

G(t-I)  r.  ( 1 )  d !  +  f  A(t-T)ltr  K  (I )  d't  (1.1) 

U,  J  -OO 

where  T  is  the  stress  tensor,  E  is  the  strain  tensor  and  G  and  A  are  respectively 
tiie  shear  and  volumetric  relaxation  moduli.  The  above  equation  can  be  inverted 
to  describe  the  strain  in  terms  of  the  history  of  the  stress.  These  forms  of 
the  constitutive  laws  are  valid  only  so  long  as  the  material  behaves  as  a 
linear  viscoelastic  material.  Figure  1  shows  the  typical  uniaxial  rate  of 
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i-'iqur*-  1.  Typical  lost  curves  for  snow  and  linear  viscoelastic  materials 


iK:  forma  t  ion  and  creep  response  cuivos  that  icsult  with  linear  vi  kcoi;  1  ant  i  c 
materials .  This  is  shown  in  the  case  of  solid  lines.  Actual  response  curves 
lor  snow  are  shown  with  the  broken  lines.  In  order  for  a  material  to  be  truly 
viscoelastic,  the  creep  curve  must  show  the  characteristics  of  the  solid  line. 

For  instance,  in  Figure  1(a),  the  creep  curve  consists  of  two  carts.  The  first 
part  represents  the  strains  which  develop  when  the  constant  axial  stress  T  is 
applied.  The  second  part  represents  the  recovery  which  occurs  when  the  load 
i.s  released.  A  viscoelastic  rebound  which  equals  the  initial  elastic  strain 
occurs  immediately  upon  load  application . 

Snow  does  not  show  this  true  viscoelastic  behavior.  Tin-  strongest  evidence 
of  this  is  demonstrated  in  the  creep  curve  when  the  stress  is  released.  fn  the 
case  of  snow,  the  instantaneous  elastic  response  may  be  but  a  small  fraction  of 
the  initial  clastic  strain.  Aside  from  this  case  the  crooi  and  deformation  rate, 
tests  do  not  readily  show  this  particular  behavior.  in  particular  the  deforma¬ 
tion  rate  tests  do  not  show  any  significant  differences  between  the  behavior 
of  snow  and  that:  of  a  linear  viscoelastic  mater  ial.  By  observing  the 
deformation  rate  tests,  one  may  erroneously  conclude  that  snow  is  indeed  a 
linear  viscolastio  material.  The  rate  curves  shown  in  Figure  1  consists  of 
two  parts.  The  first  part  of  the  curve  shows  the  buildup  of  stress  as  the 
material  is  subjected  to  a  constant  uniaxial  strain  rate,  whereas  the  second 
part  ijives  the  relaxation  of  stress  after  the  strain  rate  is  stopped  and  main¬ 
tained  at  a  constant  value.  Some  differences  do  emerge  in  the  relaxation 
curves,  but  generally  it  is  not  very  significant. 

The  actual  behavior  of  snow  exhibits  some  characteristics  of  an  elastic- 
vi scop  las tic  material,  since  some  of  the  initial  elastic  response  when  loaded 
is  jradually  dissipated  away.  However,  snow  does  also  exhibit  some  very  definite 
viscoelastic  characteristics.  This  makes  snow  a  difficult  material  to  model 
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with  a  single  constitutive  equation  even  in  the  small  strain-low  strain-rate 
region.  As  a  matter  oE  pragmatism,  most  past  investigators  have  chosen  to  model 
snow  with  a  linear  viscoelastic  representation  similar  to  that  of  Eg.  (1). 

For  most  applications,  such  representations  have  been  found  to  be  acceptable. 

The  above  is  seen  to  lx:  particularly  true  when  one  considers  other  problems 
with  characterizing  snow  properties.  Snow  is  a  highly  variable  material.  Its 
properties  change  drastically  with  temperature,  density,  and  the  effects  of 
temperature  gradient  metamorphism.  Such  properties  as  fracture  stress  and 
compliance  can  go  through  an  order  of  magnitude  variation,  and  correctly  defining 
these  properties  is  quite  difficult. 

*  i  - 1  ,  •  — 

When  finite  strains  or  large  strait,  rates  j  E  j  la  ,  t  !•;  in  ) 

arc  involved,  snow  exhibits  definite  nonlinear  properties.  Sul  in  (peg),  Brown 
••t  a  1  (1)73,  1  )74,  1-177)  ,  Desrues  and  others  (1981)  have  characterized  these 
properties,  balm  considered  primarily  the  res;  onse  of  snow  to  uniaxial  deformation} 
end  subsequently  developed  a  nonlinear  viscoelastic:  model.  brown's  later  work 
was  essentially  a  thermodynamic  formulation  in  which  compression,  tension  and 
shear  test  data  was  used  to  develop:  a  three  dimensional  constitutive  equation 
for  snow.  This  formulation  was  shown  to  be  quite  accurate  for  a  wide  range  of 
deformation  conditions.  However,  the  constitutive  equation  was  complicated  and 


■  umber soim ■ .  Therefore  the  applicability  to  engineering  problems  is  limited. 

[V’srues  and  others  have  recently  formulated  an  incremental  constitutive  equation, 
and  this  work  should  show  good  promise  for  solving  some  important  engineering 
;  roblcms . 

The  above  nonlinear  formulations  are  still  quite  limited  in  the  range  of  strain 

rates  for  which  they  are  valid.  All  of  the  theories  incur  rather  large  errors  for 

-4  -1 

strain  rates  above  10  s  .  Under  tension  and  shear  those  strain  rates  cause 


.  i  i- •<  /; i  I'-VJ 


.***'  *  — ,  -  ■ 


f racturi'  (Brown,  1974,  1977),  but  under  some  compressive  deformations,  sustained 
loading  can  be  applied  to  snow  without  causing  fracture. 

There  are  a  large  number  of  problems  in  which  the  above  discussed  constitutive 
equations  do  not  apply.  This  includes  deformations  which  produce  extremely  high 
strain  rates  (  j  E  j  10  ^s  ^)  and  large  strains  (  1 1-1  [  >  0.2).  Problems  involving 
shockwaves,  vehicle  travel,  projectile  impact  or  penetration,  among  others  are 
examples  of  such  situations.  For  instance  vehicles  traveling  over  a  low  density 
snowpaek  can  produce  strain  rates  in  excess  of  10  s  ^  and  increases  in  density 
of  loot.  Harrison  (1081)  reports  on  vehicle  mobility  studies  of  tracked  vehicles 
in  snowpaek.  In  instances  where  little  vehicle  slippage  occurs,  the  deformation 
of  the  snow  under  the  track  is  largely  volumetric,  and  this  problem  can  be 
analyzed  witli  a  relatively  simple  volumetric  constitutive  law  (Brown  1979, 

1981a,  1981b).  In  cases  where  the  vehicle  incurs  significant  slippage,  then 
a  multiaxial  constitutive  law  is  needed. 

Explosives  are  used  to  .initiate  avalanches  in  recreational  areas  and 


along  highways  exposed  to  steep  open  slopes.  These  explosives  produce 

intense,  high  frequency  shockwaves  which  propagate  into  the  snowpaek. 

Depending  on  the  explosive  speed,  snowpaek  properties  and  placement,  the 

5  -] 

shockwave  spreads  and  attenuates.  Strain  rates  on  the  order  of  10  s 

and  volumetric  strains  as  large  as  200%  can  be  induced.  Tn  order  to  analyze 

these  problems,  very  specialized  constitutive  relations  are  required. 

The  military  as  also  concerned  about  several  problems  concerning  impact 
loading  of  snow.  These  include  fuse  action,  ballistic  penetration,  and 
effectiveness  of  artillery  explosives  in  snow.  For  instance,  artillery 
fuses  which  work  well  in  mud  fail  to  detonate  properly  in  snow.  Also  rifle 
fire  has  been  shown  to  penetrate  much  less  deeply  into  snowpaek  then  origin- 


--*•  -v  a 


ally  thought.  The  unique  properties  of  snow  are  responsible  for  these 
surprising  results. 

Prior  to  the  results  reported  here  and  the  papers  published  by  the 
author  (Brown  197y-81),  many  of  these  problems  could  not  be  analyzed 
analytically.  In  order  to  do  so,  appropriate  constitutive  relations  arc- 
needed.  This  report  brings  together  tiie  work  of  the  author  on  the  high  rate- 
volumetric  properties  of  snow.  Also  included  are  applications  to  several 
problems.  Finally  recommendations  to  further  work  are  also  given. 

The  work  reported  here  concerns  only  the  high  rate  volumetric  properties 
of  snow.  Deviuforio  properties  are  not  considered  lie re .  The  reason  for  this 
is  that  very  little  is  known  about  the  high  rate  properties  of  snow,  and  a 
good  starting  point  would  be  volumetric  properties.  In  addition  some  important 
i  roblems  can  bo  studied  if  only  the  volumetric  properties  are  known.  Once 
the  volumetric  properties  are  known,  it  is  felt  that  these  results  can  be 
used  to  help  evaluate  the  do via tor ic  properties.  A  complete  multiaxial 
constitutive  equation  is  necessarily  qute  complicated,  and  it  would  appear 
to  be  best  to  approch  this  development  in  a  systematic  manner. 
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I  T  .  HIGH  KATE  VOLUMETRIC  PKUPbKTbKS  OK  NOW 

I’ract.ially  al  l  of  the  previous  work  <:j  ted  in  i: I •  int  roduct  ion  us  •>' 
purely  phenomona  logical  methods  of  characterizing  snow.  Thi  s  method  use.; 
creep  data  or  deformation  rate  data  to  determine  material  co*-f f icie»t.;  in 
the  constitutive  equation  chosen  to  represent  t  lie  mater  in  1  .  Tie-  form  of  tii* 
'•'institutive  equation  is  chosen  after  the  test,  data  is  observed  to  determine 
if  t  he  material  is  elastic,  v.i  :;c;oe  1  .c;  t  i  <: ,  linear,  nonlinear,  etc.  In  fnis 
:u  ;  roach,  litt  le  or  no  concern  is  cjiven  to  the  mic:  "Structural  de format  ion 
processes  and  their  effect  on  the  properties  of  Uv  material. 

one  reason  that,  the  phonontona  logical  approach  is  used  so  much  is  that 
it.  is  a  relatively  straight  forward  method.  Curve  fitting  techniques  such 
as  the  linear  and  nonlinear  least  squares  methods  are  commonly  used  in 
con  junction  with  statistical  analysis  to  determine  just  how  well  a  constitutive 
equation  works.  This  approach  is  very  straight  forward  for  elastic  and  linear 
viscoelastic  tna tor i a  1 s ,  but  for  elastic— vi scoplastic  materials  this  approach 
can  In  ■come  much  more  difficult,  since  very  intricate  mathematical  models 
(frown,  i  >77)  must  be  used. 

Mere  recently  Brown  (l'J/'J,  J‘>80)  ,  and  St.  Lawrence  (1977) 

,idV"  utilized  nu  crodynamicu  l  models  to  formulate  constitutive  theories  for 
snow.  This  api  roach  considers  deformation  mechanisms  at  the  crystalline 
1  •  v>  1  which  ate  rrq  ortant  in  determining  the  macroscopic  behavior.  This 
im-t  hod  of  studying  the  material  is  justifiable  in  the  sense  that  the  be¬ 
havior  of  tin*  parent  material  (ice)  and  the  dt» format  ion  mechanisms  at  the 
granular  level  must  be  responsible  for  the  properties  of  the  granular  material. 
Also  this  method  results  with  constitutive  laws  which  give  much  insight  into 
what  mi crostructura 1  changes  occur  during  deformation.  Mi crodynatnical 
approaches  would  consequently  appear  to  be  preferable  to  phenomena  log i ca 1 


formulations.  How  -ver  ,  the  mi  ci  odynumi  cu  1  approach  usually  yi'-lds  '-o:lsti- 
tutive  relations  which  are  matlnnii.it  teal  I  y  com|  1  i  '-.it  •  <].  in  fact  t  h< 
methods  have  to  date  not  found  a  lot  ul  eti'ji  noer  inq  use  tor  this  v-ry 
reason,  .tonkins  (1*»H1  ),  and  Br  own  Qy80  '  are  .  -xatni  I  os  of  r*  -cent  work  fen 
qrunul ar  material.  ;  such  as  sand  and  qr«i  1  n.  Tin.  next  decade  should  •  much 
progress  in  this  area. 

Presented  liere  are  two  formulations  devej o| ed  by  tin.  nuuini:.  '!h--  first 
method  is  applicable  to  medium- to-hi'jh  density  snow,  and  the  sc  cos  id  is  valid 
for  lower  density  snow.  Both  formulations  use  the  mi-.:rodynamicul  ;  roach 
and  hence  are  based  on  the  material  properties  id."  ice  and  the  -jr  ,»n*jiar  tux'  uri 
of  tire  porous  material,  snow. 

XI  A.  .  A  MATtKIAI.  RiiPKES  ROTATION  01'  I  CD 

Previous  studies  by  Dillon  and  Anders  land  (l;i.7)  and  Hawkes  and  Mel  lor 
(  1  ■ ;  7 1: )  have  indicated  that  po.l  ycrystal  1  r  lie  ice  under  uniaxial  tension  and 
compression  ixassesses  a  ra to -dependent  yield  stress.  Additional  experimental 

work  by  Haynes  (private  communication)  has  extended  the  previous  work  to 

_ 1 

strain  rates  in  excess;  of  1 .Os' .  In  tension,  ice  uenerally  fractures  at. 
hiifh  strain  rates,  and  this  critic:.!]  stress  is  not  st  rurally  rat-e  dependent. . 
However,  under  compression,  the  critical  stress  shows  a  definite  rate 
dependency,  even  at:  the  rates  studied,  by  Haynes .  .  Tile  present  study  primarily 
concerns  tin.*  compressive  properties  of  ice  for  intermediatc-to-hiqh  d*.  fonii.it  ion 


Based  on  tin  experimental  results  indicated  in  Fiyure  J  (Dillon  and 
Andersland  (lyf.7]  ;md  Haynes  ),  the  followinq  constitutive  relat  ion  is 
assumed  for  polycrystalline  ice: 
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,  and  A  are  material  constants;  S,  o,  and  !:  .uv  l  if  i  rii.  ipal. 

l  A 

•  1 1  ffcrenee  values  of  the  doviatoric  stress  1  onsor  ,  l  -v  iatori  c  st  i  ain  fen r.nr  , 
and  Joviatonc  rat'  of  deformation  tensor .  Similar  definitions  hold  for  D 
and  e.  The  matrix  material,  ice,  is  assumed  to  be  incompressible,  so  that 
the  de-viator  i values  of  the  strain  tensor  and  the  deformation  rate  tensors 
equal  the  total  strain  and  deformation  rate  tensors. 

The  above  equations  are  essentially  those  of  an  elastic-viscoplastic 
material ,  so  that  the  material  behaves  elastically  until  a  rate-dependent, 
yield  condition  is  readied,  whereupon  plastic  deformation  ensues.  These 
equations  are  conn  what  different,  from  <11  on's  flow  law, 

( 1 1  .  A  .  4 ) 

w : i . •  i e  and  ,  .us  ,  respectively,  the  shear  stress  and  shear  strain  rate,  <7 
a  ma'-e  ini  . instant  and  n  is  a  scalar.  liquation  2-3  were  chosen,  sins 

•  /  a,. -quit'-  a  straight-  lino  form  on  figure  1  and  fit  the  data  better  than 

•  '.lea's  iaw.  -ILeu's  lav.  does  work  well  at  lower  strain  rates,  however,  and  in 
many  asc  , mat  liemat  ical  ly  more  tractable  than  the  form  of  eqs  .  2-  J . 

insider  now  a  compression  test  on  polycrystal ! ine  ice.  The  only 

nonzero  stress  component  is  ,  where  the  x-diroction  is  the  axial  direction 

xx 

■  :  (;.•  compression  test.  The  doviatoric  stress  sensor  is 


,  ,.*n 

1  a  Cf 


(t  r  1 


(II.A.r.) 


win  re  tr  rinf  lies  tin.  trace  of  the  tensor,  1  is  the  identity  tensor,  and  •  : 
is  tlie  t.'uchy  stress  tensor.  S  therefore  has  the  <  um]x>nents 


10 


11 


-■  s  M.nn  ( '■  u  )  i ;  •  .a  .  u } 

XX  0  ,!  XX 

by  .id  jilMl  inu  '■<[.  ].  I  to  I  It.  t_  i  i<  •  com;  t  •  ■  .  i  on  d.it.-i  ;■•:!.  .  u  J  .  j  t; .! 

A  ;  lx  lit  N  IP 

1  .  lfoxlO  1:  “  •  !  T  .  ,\ .  i .  j 

A  -  i.  lxio’ 

As  can  be  seen  in  Figure  2,  this  equation  re;. resents  the  n-su  1 1  s  quit  •  • 

w.  11.  There  is  an  appreciable  amount  of  scatter  indicated  by  the  v<  rt.i  m I 
lines.  This  is  to  be  expected,  as  there  may  be  a  qood  degree  of  variability 
in  freezing  the  ice,  in  forming  the  specimens,  and  in  testing.  However,  the 
average  results  appear  to  be  well,  approximated  by  the  constitutive  law. 

The  constitutive  law  does  not  involve  a  work -hardening  term.  Previous 
expel  imental  studies  on  ice  by  Dillon  anil  Anders  land  ( 3.062)  ,  among  others, 
have  not  indicated,  strong  vork-hardeni ng *  characteristics  tor  ice.  However, 
there  are  no  published  data  on  ice  involving  large  strains  in  excess  of 
several  hundred  percent.  Therefore,  it  is  quite  possible  that  at  large 
strains  ice  does  exhibit  work-hardening  properties.  Such  strains  certainly 
do  occur  dui i ng  finite  compaction  of  snow,  particularly  in  the  very  critical 
grain  bonus  where  massive  localized  deformations  take  place.  Thus,  the 
constitutive  equation  finally  formulated  for  snow  may  have  to  bo  empirically 
adjusted  with  a  work -harden j ng  term. 

This  problem  has  also  been  encountered  in  the  field  of  powder  metallurgy 
where  a  porous  metal  witli  an  elastic,  perfectly  plastic  matrix  material 
shows  considerable  work-hardening  characteristics,  yuite  possibly  the  large 
strains  occurring  in  high  stress  regions  such  as  grain  bonds  are  of  such  an 
order  of  magnitude  that  the  matrix  material  does  work  harden.  Also,  quite 

* A  Work-hardening  material  is  one  that  stiffens  under  plastic  deformation. 

U 


possibly  Mlt:  actual  di.-sei  i  I  t  oil  ol  l  ;i<  let  anal  ion  :■  I<  .  !■:  i  :!■  '  >ri;  .  /it 

any  rati-,  the-  work -hardening  hat  j/Ltr  j  st  ol  porou  ;  mu  t  ■  t  •  i!.:  ,, 

;  'vd'  Il.1  aluminum  arc  not:  predicted  by  i*x|  tr  i  nunt  s  t  nut.  have  n  r  in  or  the 

solid  Matrix  material. 

I  1.1..  I 13  V  bhc/1  WENT  OK  THU  VOUIMETK 1C_ -JUNUri Tl !T1 VE_ U»W  KOH  K*K  i'U.hAl’EE  Mt  bi: 
Consider  now  the  deformation  of  a  thi  ck-wa  I  !•  d  hollow  r.jivt'e  of  at.  i  i , 

■  ••.mprossible  viscuj  last  ic  material  with  the  constitutive  relation  qiven  bv 

■  lr.A.1  -  II. A.  3.  The  ueometry  is  jiven  in  Figure  3.  The  initial 

internal  and  external  radii  of  the  sphere,  a  arid  b  ,  art  choser.  so  that  the 

o  n 

correct.  material  porositv  results.  1  ri  this  paper,  the  density  ratio  is 
defined  as 


-  •  V  ■ 


rutrjrat ion  yields 


'  -Htt  ) 


m  .  b  .  4 ; 


whore  K  i  s  an  unknown  function  of  time.  This  ion  can  different  i..;'"d 

twice.  Th i s  ):•  ■  --u  1  *  a  i  n 


for  tile  radial  acceleration.  The  acceleration  can  be  expressed  in  terms  of 


■deration  potential  i-(r,t): 


( I  T  .  B  .  6 ) 


by  do  ina  this,  the  fol  lowimj  can  be  arrived  at  for  all  a-'  r  <  b: 


B  Fi 

.  d'  +  3r ' 
lBr 


(II. B  .  7 ) 


The  above  results  art:  strictly  of  a  ki.ncmati.cal  nature  and  depend  only 
on  the  constraint  of  material  incompressibility.  In  addition,  material 
i  :i'  ’ompross  ibi  1  i  ty  may  be  used  to  arrive  at  the  fol.iowiny  relations: 


b  .  3  J  J, 

‘o  “  V(W 

(b  -a  ) 


a'!('(-l)  ,  a^t 

». _ .3  _  0 _ 

<ao"U  '  "  ('Vn 


•'u 

-vir 


(II .B.Ba) 


(II . B.Sb) 


(II . B.flo) 


As  the  external  pressure  1’  increases,  the  deformation  proceeds  in  three 
distinct  phases: 

(1)  an  initial  purely  elastic  phase 

(2)  an  elastic-plastic  phase  with  an  elastic/plastic  interface 
at  r  =  c,  where  a  <  c  <  b,  and 


(3)  A  fully  plastic  phase .  During  the  first  two  phases,  the  strains 
are  assumed  to  be  small,  but  dur  inn  the-  fully  plastic  phase,  large  rains 
can  be  incurred.  Kuch  phase  is  now  considered  separately, 
i  ul  ly  F’hasi  : 

In  this  phast  ,  the  strains  are  a  ,sutr,  d  to  be  ini  i  ni  tesimul  .  "'h<>r.  lore 

'he  three  strain  component;;  are 
■u 

*  ,  ,  -•  >  ,,  -  u/  r  ill.  !•, . 

r  !  r 

where  u  is  the  radial  displacement;  and  if  u  is  small  compared  with  r ,  u  - 
r-r  may  be  approximated  by  the  expression 


u 


HJ_t_) 


(II. B . 10) 


For  small  strains,  r  differs  by  a  small  amount  from  r;  therefore,  in  this 
section  the  distinction  will  be  dropped.  The  strains  become 


2  b 


( 1 1 .  B  .  1 1 ) 
(II.B. 12) 


Since  the  material  is  as:  jmed  to  be  linearly  elastic,  the  constitutive 
i.  .;uation  acquire.1;  the  form 

(II.B.llal 
(II. B. 13b) 
( 1 1  .  B  .  1 3c ) 


1  -  2Ge  =  4d  B/(3r  ) 

r  r 

20,.  =  -2(111/ (Jr*) 


S,  --  2i'ie.  -  -2GB/ (3r  ) 

'■P  9 

where  S  ,  S,.,  S.  and  e  ,  c  e.  are,  respectively,  the  deviatoric  stress  and 
r  b  <p  r  0  (J) 

deviatoric  strain  components.  Since  the  material  is  incompressible,  the 
volumetric  strain  is  zero,  and  the  hydrostatic  pressure  p  cannot  be  evaluated 


with  tlie  constitutive  equation  alone. 


-NrtJiKj... 


The  radial  equation  ol  motion  i : 


: - t-  (■’  -0,1  "  •)  I 

■  )'  r  r  m 


i ioh  suv  tie  nut  in  the  modified  ton 


![!.!'.]  -1 ) 


■  i  m  r 


( i  r .  i-. . !  h ) 


or  whit;:;  •  bou;  .’a /  conditions  are: 


•id.)  r  --  h. 


(r  t  ir>) 
(II -B- 17) 


should  be  noted  that  pore  pressure  is  assumed  negligible.  Under  extremely 
large  rates  of  loading  or  near  low  <1.  -nsity  rat  ios  (say  !«'<<  1.2)  ,  tho  wisdom 


a  .Mira;  tjvJi  may  1 e  que.stionjMc.  Tin:, 
.  Integration  el  eg.  Is  result:;  in: 


-ill  ,1.)  .  ,  (r , t)  ih  ( 1 ) 

in 

•vi.ej-e  in*)  is  a  time  f  unct  ion .  Since 


<n  2  J  be  discussed  in  more  detail 


(IT .B . 18 


+  I,  +  4GB  (3r  ). 


■irvi  t.lie  boundary  conditions  gives  the  solution 


(1I.B.19) 


:t. ) 


4GB,  1  1 

+  (  -  -  -V,- 


I  ,  7  > 


(T  1  .B. 20) 


'  can  then  be  expressed  directly  in  terms  of  l'(t)  by  making  appropriate  use 
of  pcs .  7  and  8.  These,  combined  with  eq.  20,  yield 


where 


2  .  •  4G<a0-u) 

1’  (t)  -  1  QUx,u,n)  +  ^ 

3a((X-l) 
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(  1 1 .  R  .  2 1 ) 


(I  l  .  H  . 2 2 ) 


Fount ion  22  describe;-  tin:  material  volume  t.r  If;  r< -s:  ions':  t.o  a  iiydrost.it  Lc 
;  ai  loading.  ! '  i  i « ■  term  I’ ( l  >  i  r<  ^r,1  a  t: ; actual  j  »r«  o  in  the 

ii.ati'ix  material .  Cat  rol  1  and  Uo.lt  ( 1  ‘>72)  ha-/"  shown  that  the  average 
!  rtssuie  in  tin.  porous  material  p(t)  i«  approximately 


P  =  I’/  « 


( 1 1  . 1'. .  2 1 ) 


At  lony  at  tin:  material  does  not  yield,  oq.  11  doscr ibes  the  time— 
depend. ait  response  of  the  material. 

This  ;u.it  ion  is  the  same  as  Carrol  l  and  Holt’s  (  llJ72)  result,  in  ir 
;t  udy  f  1  jorous  a.luminam. 

! 1  i  a  .  t  it:-!  ■ !  a  . ;  t  i  r  i  has.  : 

equation  l"  indioat.es  that  the  maximum  stress,  occurs  at  tin:  inner 
i.i.iins.  V  i  el.il  no  then  ii.it  Late:.;  I.liere;  and,  as.  tie  pressure  I'  { t  )  rent,  i  nues 
it.  i  net  • -as..  ,  a  yi.-ld  surface,  r  ■•(  t  )  ,  pro;  -agates  radially  outward.  The 
yi.-ld  cot.di  t  Lon  gi  von  i>y  eg.  (II. A. 3)  is  reached  at  a  critical  pressur.  }■ 
when  the  principal  deviatoric  stress  difference,  S  -S^ ,  reaches  the  critical 
value  Y,  which  itself  is  rate  dependent. 

Assume  now  1*  •  !•’  and  that,  yielding  tias  propay  a  tod  out  to  a  radial, 
distance  e,  a  ' c  b.  For  r  •  c,  an  inelastic  stress  state  exists,  and  for 
r>c,  the  material  is  elastic.  When  analyzing  the  outer  elastic  zone,  the 


following  boundary  conditions  apply: 


o  =  -r  (t)  ,  r  =  b 

r 


(IT.B.24) 
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i 


p  ( t- ) ,  1 
c 


;b 


(  T  1  .  H 
(  I  1  .  I; 


At  the  interface  ui  the  clastic  zone  and  elastic  zoic  ,  tin 
satisfy  t.he  yield  condition.  In  the  yielded  zone,  •.<;  1!., 


st  fee s  mu:  t 


\ .  5  lj  i  v<  s 


In  A U) 

‘  0 


!  i  ! . 


'.ciLoh  remains  valid  as  long  as  the  loading  is  mono  tonic. 

consider  first  the  yielded  zone,  whore  the  constitutive  lav.-  in  defined 
lii  terms  of  the  deformation  irate1  tensor.  The  principal  rate  com:  unents  arc 
:  t  and  .  arc-  t o , 


I)  ■■  -r 

i  r  r 


(  I  I  .  b  .  2*) 


1 i 


( I  L  .  H  .  2 '  i ) 


.  ---  '  ,/r  (l  2i  ) 

r1  V  * 

shear  in.:  components  are-  zero.  since  lie-  deformation  is  viscous,  a 
l.ai's'ian  description  can  lie  used  i:eic.  1  ncoi:i|  .-m  •::::  i ! i  I  l  t.y  require::  that,  the 
>ii  v. -r  nee  of  the  velocity  Vector  vanish;  j.e., 


— (v  r  )  =  0. 
r  r 


( 1 1  .  B .  U ) 


T!  lie  yields 


K  ( O 


r  2. 

r 


(U.B.  32) 


where  F  is  an  undeterrai nod  time  function.  Differentiating  the  above  and 
equating  to  the  accelerat. ion  potential  results  in 


F  1  I-'_ 

--  •>  A 


{  r  I  .  B  .  .3  3 ) 


Actually  !•’ (t)  can  be  related  direct  I  y  to  h  { I  )  •  1-  - 1  l  u*  ••  1  in  cj.  j  .  7  ,  which  mu 
!"■  Join  •  by  differentiatin';  B  f  o  obtain 


IM  t  )  - IK  ( t ) 


(  .  i  .  P.  J4) 


chic!)  ultouJ.il  be  expected,  since  F  and  B  have  lieeii  Lased  on  the  ju.-ely 
k  i  itemat  Leal  constraint,  of  material  incoinpressi  bi  1  i  ty . 

Utilizing  the  incom;  lessibi  1  ity  constraint:::  in  •  H  -jives 


F(U  -  V* 


(I  I .  B.  J!>) 


3lvu 


lev:  idr-r  t . :  t  -  -  racii  n  equation  of  motion.  in  frris  of  too  <leviat.fr  : 

t.  »  .-SS  coin:  eneh!  J  , 


•r  tn  <r 


11I.E. 


rincital  -viator  ic  cii  format  i  * .  rat..-  <ti  *  !’•  c-nt 


-•j  -  if  •  r 


.  !  :-B.  17; 


l  !,-  .lit  f-  r'-nci  --S  'an  easily  !»■  m  -  n  to  iiav--  the  vain- 


o<  ••  i  '.tr  in 


t  in  ■  an  t  ’  r  mo'  ion  '.jives 


■  <  t  )  •=  ”-S  +2ii.  In  (r/b) 


'  1  n  (  r/r )  [  -  •  1  n  ( r/:',)  J 


{  1  I  .1:.  IS') 


,  ^  k  (t) 


--  [  -  3AF ! 


(11 .B. 17) 


m<i  k(t)  is  a  tint*  integration  function  resulting  from  solution  of  the 
■quatioii  of  motion.  if  we  denote  as  1’  the  value  of  the  radial  stress 
:t  tb<  interface  r  --  c,  appl Lea t ion  of  the  stress  boundary  to  eq.  3B 


P  :  2(S  -C)ln(c/a)  +  3C1  n  (c/a)  In 


where  and  ^  denote  the  value  of  tji  at  r  =  a  and  c,  respectively.  T 

must:  th''ii  be  -.•quoted  to  the  radial  '-i  u  :!s  obtained  from  tin;  sol:::  ion  o 
-  1  as  t:  i  e  zone  C"  r  •  1  - . 

In  t:hv  elastic  zone,  a  :  voc  dur>  ■  i  d- -nt  i  e  i  a  1  to  that  follow'd  ii.  t. 
first.  ;  ha  so  <-t  th-  dot"  irinati-a:  c-iii  i-i  followed  to  arrive  at  !::••  !  - ,  1  low 
soii.d  ion: 

!  i  -  eh  ( c  ~u  ' 

1  n  . ;  :  j 

i  .:,;n  b.  e  i  i  in :  noted  from  eqs.  -S'  ,j,.j  4r  : 

t  ' 

■(f)  N  d.  (  "-b_  ')  -  1  n  (  "-)  (2S.- U'lna  1 

a  o  -If. 

:  '  i . 

.-  -:u;i  : -ui-s- -queut  1  y  e  1  i i  nat  ed  by  means  of  t-qs  .  ,  17,  41 

.a:-:  so  that,  -i  re  la  t  iotmh  i between  i  and  !  c<lli  bo  established, 

.'in.--,  i>,  -.lil  !..-■  shown  la  tor,  this  i  n  termed  ia  t  e  phase  is  riot,  of  much 

;  at  .  t.nij-  final  equation  'will  not  bo  developed  here, 

ru  •  1"  l  lriiit.  i-:  phase: 

ere-  reaches  b,  the  sphere  becomes  fully  viscous,  and  the  full, 
u  r-  -.: , i : .  !«•  considered  with  one  constitutive  equation.  The  equation 
of  i-ui,  once  the  appropriate  substitutions  are  made,  once  aqain  vie 
.-  : .  :K,  .-xc-.-j-f  that  this  time  tin.  Iioundary  conditions  qivc-n  by  eqs.  if- 

17  -q  ;  1 y ,  Aj  i  licat ion  of  the  boundary  conditions  results  in: 


]  /l 

Pit)  -  In  If—  — )  )  2  (S  -C)  t-Jc'ln 

(-1  0 


■;i 'I'f’rMiAt#!  ■ 


As  indicated  by  4.;,  tile  i  ressure  1 e:,j  a  m:  ■  of  t  Ir-  malarial  c  ,i.  •  l  : '  . 

two  i'.nts:  (a!  a  ou.i :  i ;  t  at  i  c  part  an-l  (1.)  a  dynamic  part  as  t..  :  r  •:.!  ■  d 

ov  tin':  ,ta  .  li-rat- 1 ,  in  ; ...  _■>  t  ■  ntia  1  t<  -rut  .  'C~.  )  /•  i . 

in  a 

The  a  ci  ■  1  ■  'fit  ion  r-  rm  cat;  ensile  shown  to  iji 


,  )  -  - ,  ( 

a  Jr-'-,- it  b 


(b  ‘‘-a  ') 


i  -  J  ) 
c 


...  .  .tin:  b  at  i  •  x;  •  a  i  i .  i  •  •  in  u  r;\;  ■ , •  by  •  ■  -i . 

ar  •.•!«  ra*  nil  t  •  •  rw  becunvc.  niuni  f  i •  ■  :.!•,•  at:  vry  at  b  aa  r  c  • 
w  tch  inv..  lv"  si  i  f  i..ii;t  inertial  •  ff.  .-Is,  includ-isr  shook  v.av  s ,  i  :n|  ret  , 
■  t  ■  • .  :•  ♦•••sent  !y  v.1.  .if.  •oti;  i-rif  <t  vci  t:  ti  qaa.;:  -s'  at  ir  1  ...!  s  :  ‘  .  t  .or...  that 
t  ■  . i . a ■  1  ra'  ion  t  ••rtn.-t  may  be  ta  il.  'of  ■  i. 

As  i  r..li  c.itt.-cj  section  11. A,  modi  f  ic-'it.ioii  o!  Vito  ..<;ns  t:  i  :a  t  i  v  • 

law  is  nocdotl  to  r.  f  1  .-•vt  wot  k-'nat  deninq  of  foots  and  to  provid-  s.on.e  .-or  - 
-t  !■  ti  to)  sotiii-  tin-  : . ;  t:.pl  i  !  y  i  nq  a :  :  a  tin;  / 1  ions.  us-d.  l.y  idoinq  >  j---  1 1  .. 
to  •rK-bartloiiiriff  t  erm,jex|-( -'*>•«/•<  )  ,  and  by  makinq  ns.,  of  24,  ..q  'J  dun  i  • 


put  in  tits  followiny  form: 

■  ti  ln( ;:\!l 


i  A)  i 

:  in  i  •  i  .)  "xi  (-;  ,/  .  .  (u.i' 

‘  t  ■  - 1  i  > 


r>  y  and  J  .  a  r  ••  mat  r  iui  ooot'i  icL.-nts.  T'u  i  work -hard' -t.i  tiu  is  ■!"  a 

:. atari.  . 'onsistent  'with  t  h>  •  term  used  by  St  .  l.awreiico  and  brad  ley  (b'”;4)  , 
who  not’d  that  oilman  (!"(  ’)  used  a  similar  term  tot  lithium  fluoride 
irvstals  to  account  for  vari.at  .otit  in  the  mobile  dislocation  Velocity.  It. 
porous  snow,  work-hnrdi.  n  j  nq  may  lit-  due  t  o  <  xo-ssi  V"  amount  a  of  iutcriir.iiiul.il' 
motion  that  would  not  be  present  in  solid  i  olyt  rystn  1  1  i  nu  ice.  St  .  Lawrence 
and  fu.idjey  (  ly74)  addressed  this  question  in  a  mote  thoiouqh  manner.  In 


22 


'•  »i  ■■!+•&%,  id- 


add i  I. i  on  ,  tin1  mod'  •  i  used  :  j  i  ■  a  •  -  is  t.ifii.i  id'll  ii.  iii'  'iii'-  t  r .  i  *  i  os. 

of  till'  [  ores  j  .ISSIlinod,  SO  tll.lt  t.!<:  oil  .-t  of  OH"  lot"  ol,  dill  » t  v..i  . 
“t'jl.  t  d.. 

S 1  ;;i:  1  i  t  l .  i  iu.it.  :  i  ■!:  : 

if  Y  is  coi.sidornblY  stitct  L I  ".r  t':,.i:i  d,  it  o  Hi  i  ■  shewn  Vs. a'  ‘ 
i  ii.-’trt  f.l  ( and  t  li ■  •  i ■•■t<>t  •-  t  ,)wi  1  I  •  'ju  :  wan !  I  'im  it~  t  !.«•  i  i  •  t  o 

:  iu''fs  of  tlu  •1i-f('*-::.dt-  iuti.  If  «wh  is  td-s  .;.a  t.ii- •  w  .....  s  ot 

:i  :  r::i,-»t  io:.  cm.  :  ••  o"  Vy  lx-j  r>.  J--:.*  in  i  fi.o  m.ilitL.!  ,.s  i  r  i  :  .  :- 


■  ■  d-;;  Ar  i.  idv  ;  ;  x;  idi  ii-’riN’iVf;.  data  kok  i-om:  I'obL/U'.Sf  mme;. 
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The  density  ratio,  however,  is 


a  -  a0tl-(X/h0)t] 


(II. C. 5) 


and  its  rate  is 

a 


(II. C. 6) 


The  constant  rate  tests  do  give  a  constant  density  ratio  rate.  Note, 

however,  that  the  rate  of  change  of  density  is  not  constant.  For  this  set 

of  experiments,  a  was  readily  calcualted  by  eg  6  and  eg  II.B.47  was  then 
utilized  to  find  the  variation  of  p  and  a  for  a  given  rate.  Figure  5  shows 

the  comparison  of  theory  with  experiment  for  snow  at  -10#C  and  for  three 

i 

initial  densities.  The  data  acguired  by  Abele  and  Gow  (1975,  1976)  included 

a  variety  of  strain  rates,  since  different  specimen  sizes  and  crosshead 

speeds  were  utilized.  The  data  shown  here  reflect  rates  in  the  neighborhood 

•  <v, 

of  a’v  10s  . 


The  data  available  from  Abele  and  Gow  (1975,  1976)  measured  only  O^, 
the  major  principal  stress,  whereas  the  hydrostatic  pressure  p  =  -1/3 (o^+o^  +  ° ^ 
Has  needed.  Therefore,  a  series  of  experiments  was  run  to  measure  lateral 
stress  as  well  as  axial  stress  so  that  Abele 's  and  Gow's  data  could  be 
adjusted  to  reflect  p  rather  than  It  was  found  that 

0.6  0X<  p<0.9801  (II. C. 7) 

for  all  the  experiments  run,  thereby  fairly  well  bracketing  the  data  report 
by  Abele  and  Gow  (1975,  1976). 

One  interesting  result  was  observed  in  the  testing  program.  In  tests 
involving  "old”  snow,  the  lateral  stress  was  found  to  be  about  90%  of  the 
axial  stress,  but  for  tests  of  "new"  snow,  the  lateral  stress  was  only  30- 
40%  of  the  axial  stress.  In  each  case,  unbonded  snow  (sifted  within  two 
hours  of  the  time  of  testing)  was  used,  so  that  the  difference  between  the 
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Figure  5.  Comparison  of  pore  collapse  theory  with  experimental  data. 
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lateral  stress  in  the  old  snow  and  that  in  the  new  snow  must  have  been  due 
to  differences  in  crystal  structure,  primarily  crystal  shape.  The  old  snow 
had  undergone  equi temperature  metamorphism,  and  the  crystals  had  a  spherical 
shape.  The  grains,  therefore,  were  capable  of  rolling  and  sliding  relative 
to  each  other,  thereby  accormnodating  a  lateral  motion  to  produce  a  large 
lateral  stress.  New  snow,  with  its  complicated  grain  structure,  would  not 
allow  this,  thereby  producing  a  smaller  lateral  stress. 

As  can  be  observed  from  Figure  5,  the  comparison  between  theory  and 
experiment  is  quite  good,  with  the  essential  characteristics  of  the  deformation 
being  represented  by  the  theory. 

Figure  6  compares  theoretical  pressure  curves  with  data  collected  and 
summarized  by  Mellor  (1974) .  This  figure  contains  results  of  laboratory 
studies  as  well  as  field  data  relating  density  to  gravity  reasonable  compared 
with  the  data.  It  should  be  remembered  that  some  of  the  data  shown  in  this 
figure  represent  uniaxial  stress  conditions,  and  that  the  actual  hydrostatic 
pressures  are  only  one-third  the  values  shown  for  these  data.  Therefore, 
some  of  the  experimental  curves  would  move  down  vertically  relative  to  the 
theoretical  curves.  However,  since  there  is  such  a  diverse  range  of  load 
histories,  temperatures,  and  time  ranges  contained  in  Figure  6,  any  meticulous 
adjustments  would  not  necessarily  change  things  that  much.  What  can  be 
said,  though,  is  that  eq  II.B.48  appears  to  be  functionally  correct  for 
snow  with  initial  densities  exceeding  300  kg  m 

Figures  7  and  8  further  describe  some  important  properties.  The 
deformation  rate  dependency  is  illustrated  in  Figure  7.  For  snow  with  an 
initial  density  of  3S0  kg  n  \  the  stress  response  is  shown  as  a  function 
of  density-ratio  rate  at  three  different  instantaneous  densities.  As  can 
be  seen,  a  rate  dependency  does  exist,  but  the  importance  of  rate  decreases 
as  rate  increases.  At  rates  characteristic  of  stress  waves,  therefore,  one 
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Comparison  of  pore  collapse  model  with  data  supplied  by 
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Figure  7.  Variation  of  prfssure  response  with  rate  of  change  of 
density  ratio,  a. 
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night  be  able  to  assume  a  constant  yield  stress  and  achieve  a  simplified 

version  of  the  constitutive  law  given  in  eg.  II.B.48.  However,  for  lower 

-5  ,  • , 

rates,  say  in  the  range  10  <ja|<10,  the  rate  dependency  is  significant 

enough  that  such  a  simplifying  procedure  would  not  be  recommended. 

Figure  8  gives  an  indication  of  just  how  effectively  snow  can  absorb 
energy  during  compaction.  In  particular,  this  figure  shows  the  work  required 
to  compact  snow  to  a  terminal  density  of  700  kg  m  ^  for  a  range  of  initial 
densities  find  density-ratio  rates.  One  can  see  immediately  that  initial 
density  has  a  dramatic  effect  on  work  required  to  compress  snow.  Density- 
ratio  rate  is  also  significant  but  certainly  much  less  so  than  initial 
density . 

II. D.  DEVELOPMENT  OF  CONSTITUTIVE  LAW  FOR  NECK  GROWTH  MODEL 

The  deformation  of  porous  material  such  as  snow  involves  several 
different  mechanisms  which  act  simultaneously  to  varying  extents.  Three 
such  mechanisms  are;  (a)  pore  collapse,  <b)  intergranular  glide,  and  (c) 
inelastic  deformation  of  intergranular  necks.  Pore  collapse  and  neck 
deformation  both  involve  deformation  of  the  matrix  material,  whereas  inter¬ 
granular  glide  is  a  function  of  bond  strength  and  surface  frictional  properties. 
Pore  collapse  is  a  predominant  mechanism  at  higher  densities,  and  was  shown  to 
accurately  define  mechanical  behavior  in  the  previous  section. 

For  low  density  snow,  the  volumetric  deformation  is  determiend  largely 
by  the  intergranular  glide  and  inelastic  deformation  of  the  grains  and 
necks.  Under  a  state  of  hydrostatic  pressure,  one  may  expect  effects  due 
to  intergranular  glide  to  be  reduced  to  a  minimum  relative  to  inelastic 
deformation.  Due  to  the  random  nature  of  both  grain  geometry  and  grain-to- 
grain  bonding,  glide  effects  certainly  cannot  be  eliminated,  but  it  may  be 
reasonable  to  ignore  them  in  comparison  to  the  volumetric  strains  due  to 
plastic  deformation  and  flow  of  the  matrix  material. 
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In  this  section,  a  constitutive  law  for  volumetric  deformation  of  low 
density  snow  is  formulated  by  using  a  neck-growth  model.  In  order  to 
account  for  intergranular  slip  effects,  an  empirical  adjustment  of  the 
equation  is  then  made,  and  the  results  are  then  compared  to  existing  experimental 
data. 

From  Figure  9(a),  a  good  idea  of  the  grain  and  neck  geometry  can  be 
obtained.  The  grains  are  those  structural  elements  which  have  free  surfaces 
which  are  predominantly  concave  inward,  toward  the  center  of  the  grain, 
where  this  radius  of  curvature  is  denoted  by  R.  Individual  grains  are 
connected  by  necked-down  regions  whose  surface  are  concave  outward  toward 
the  air  phase  in  the  section  view  shown  in  Figure  9(a)..  This  radius  is 
denoted  by  R' .  The  bond  is  the  plane  containing  the  minimum  cross  sectional 
area  of  the  neck.  When  two  ice  grains  are  brought  into  contact,  the  length 
L  of  the  neck  begins  to  increase  from  its  initial  value  of  zero,  since 
sintering  effects  begin  to  produce  the  concave  outward  geometry  of  the 
neck.  This  is  illustrated  in  Figure  9.  In  many  cases  the  length  of  the 
neck  can  be  quite  large,  particularly  for  low  density  snow  in  which  some  of 
the  necks  have  developed  from  the  fragile  needle  structure  of  fresh  snow. 

If  one  is  considering  fracture  strength  of  snow,  the  bonds  are  of 
primary  importance,  since  these  are  the  regions  of  low  cross  sectional  area 
in  the  material.  However,  characterization  of  the  deformation  properties 
of  snow  must  consider  the  entire  neck  and  the  ice  grains,  since  they  con¬ 
tribute  to  the  total  deformation  process. 

In  developing  the  constitutive  relation,  inertial  effects  are  assumed 
to  be  negligible.  The  major  portion  of  the  deformation  procss  is  assumed 
to  take  place  in  the  necked  regions  of  the  granular  structure.  The  grain 
bodies  them  selves  eventually  begin  yielding  as  the  pressure  increases 
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a.  TYPICAL  GRAIN  and  NECK  GEOMETRY 


b.  IDEALIZED  GRAIN  GEOMETRY 


Figure  9.  Representation  of  grain  and  bond  geometry  for  neck-growth 
model  of  snow  deformation. 


during  the  deformation  process.  In  undistrubed  snow  that  has  had  several 
days  to  sinter  at  temperatures  of  about  -5#C,  the  grain  bonds  will  usually 
have  cross  sectional  radii  at  least  25%  as  large  as  the  grain  bodies  themselves. 
This  rate  of  bonding  depends  to  a  large  extent  on  the  temperature  and 
imposed  pressure.  At  temperatures  below  -  30°C,  snow  sinters  very  slowly. 

In  the  absence  of  a  temperature  gradient,  snow  at  -5°C  can  be  fully  bonded 
within  one  day  with  grain  bond  radii  on  the  order  of  10%  of  the  grain  body 
if  moderate  pressures  are  imposed  on  the  material. 

With  reference  to  Fig.  9,  the  ice  grain  can  be  divided  into  two  regions. 
Region  1  is  that  portion  of  the  grain  body  upon  which  lateral  surface 
loading  is  imposed,  whereas  Region  2  has  a  load-free  lateral  boundary. 

Region  1  does  have  some  stress  free  surface  since  the  grain  bonds  do  not 
completely  envelop  the  surface.  The  material  is  assumed  to  be  completely 
yielded,  since  primary  interest  is  in  snow  under  large  volumetric  deformations. 

The  grain  radius  is  R,  and  A  is  the  neck  radius,  the  bearing  stress 
P  (Fig.  9(b))  on  the  grain  bond  is  assume  to  be  related  to  the  macroscopically 
measured  pressure,  p,  by 

P  'V  pa(R/A) 2  (II .D.  1) 

Here  a  is  taken  to  be  the  density  ratio, 

a  =  Pm/p  (II. D. 2) 

where  p  anddp  ana,  respectively,  the  densities  of  the  matrix  material  (ice) 
in 

and  the  porous  material  (snow) . 

Consider  first  Region  2  of  the  ice  grain.  Assuming  surface  tension 
effects  are  negligible,  its  lateral  surface  is  completely  stress  free.  The 
stress  tensor,  using  a  cylindrical  coordinate  system,  has  the  nonzero 
stress  components  Ogz>  0^,  Ogg,  and  O^.  The  stress-free  boundary  con¬ 
dition  on  the  lateral  surface  can  be  used  to  obtain: 
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(II. D. 3) 
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0  3r/3z 

zz 


The  z  and  r  components  of  the  equation  of  equilibrium 


3a  /3z  +  (1/r) 3 (r  0  )  /3r  =  0  (II. D. 4) 

zz  rz 

3a  /3r  +  3a  /3z  +  (a  -oQa)/ r  =  o  (ii.d.5) 

rr  rz  rr  oo 

The  condition  of  incompressibility  requires  that 

3v  /3z  +  (1/r) 3 (rv  )/3r  =  0  II. D. 6) 

Z  XT 

where  v  ,  v  and  vQ  are  the  velocity  components  in  the  r,  z  and  0  directions, 
r  z  o  * 

Integration  of  eq  4  and  use  of  Leibnitz's  rule  and  the  mean  value  theorem 
results  in 

R  _  _ 

2irr(3a  /3z)dr  =  A  3a  /3z  +  (a  (z,t)  -a  (R,z,t))  3a  /3z  (II. D. 7) 
zz  z  zz  zz  zz  z 

where  a  is  the  average  axial  stress  on  the  cross  sectional  area  A  at  the 

ZZ  Z 

location  z  in  Region  2.  The  following  relation, 

R  _ 

2irr(3a  /3z)dr  =  A  3 a  /3z,  (II. D. 8) 

zz  z  zz 

is  approximately  correct  if  either  (C  (  z,t)  -  0  (R,z))  or  3a  /3z  is 

zz  zz  z 

sufficiently  small  to  make  their  product  small  in  comparison  with  the  other 
terms  in  eq.  7.  Assuming  this  approximation,  the  z-component  of  the  equation 
of  equilibrium  becomes,  after  application  of  the  boundary  condition. 


3c  /3z  +  (20  /R)  3r/3z  =  0  (II. D. 9) 

ZZ  zz 

Utilizing  material  incompressibility  and  the  predominantly  axial  state  of 
stress,  the  flow  law  given  in  eq.  II. A. 1  acquires,  after  some  algebra,  the 
form 

O  -  2S  +  2C  In  (  (3A/4)  3v  /3z)  '  (II.D.10) 

zz  o  z 
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Solving  for  3v  /3z  gives 
z 

3v  (z,t)/3z  =  (4/3A)  exp  (0  /2C  -  S  /C)  (II.D.ll) 

z  zz  o 

which  can  be  shown  to  satisfy  the  equilibrium  equation,  eq.  9.  The  in¬ 
compressibility  condition  can  be  used  to  show  that  the  rate  of  change  of  the 
radial  velocity  v is 

v  =  -  (r/2)  3v  /3z  (II.D.12) 

i  r  z 

Setting  r  and  v  equal,'  respectively,  to  A  and  A  in  eq  12  and  substituting 
into  eq.  11  yield 

A  =  -<2A/3A)  exp  (a  /2C  -  S  /C>  (II.D.13) 

zz  o 

Now  consider  the  deformation  in  Region  1  of  the  ice  grain.  In  this 
region,  part  of  the  lateral  surface  is  constrained  by  the  bearing  pressure 
of  adjacent  grain  bonds.  In  this  case,  at  a  neck 

°rr  =  "  P  a(R/A)2  (II .D. 14) 

whereas  a  is  zero  on  other  parts  of  the  lateral  surface  not  in  direct 
rr 

contact  with  a  neck.  The  average  stress  on  the  lateral  surface  must  be 
reduced  by  a  factor  f ,  0  <  f  <  1,  which  gives  the  fraction  of  lateral 
surface  in  contact  with  a  neck.  Therefore  the  average  value  of  0  at 
r  =  R  is 

0  =>  -  f  p  a  (R/A)2  (II.D.15) 

rr 

which  gives  the  mean  radial  component  of  stress  acting  on  the  radial  surface 
of  the  grain  in  Region  1.  The  radius  R  is  relatively  constant  with  respect 
to  z  in  this  region. 

In  what  follows,  the  superscripts  I  and  II  are  used  to  refer  to  regions 
1  and  2,  respectively.  From  eq.  1  the  axial  component  of  stress,  O  ,  in 
Region  1  is  approximately 
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These  last  two  equations  define  an  average  stress  state  in  Region  1. 
The  maximum  difference  value  S  in  Region  1  is  then 


1  l  (R/A) 2  f-lj 

22 


azz  11  [  f  -  (A/R)2] 


(II.D.17) 


Making  use  of  incompressibility,  the  maximum  difference  value  for  the  rate 
of  deformation  tensor  can  be  found  to  be 


D  =  (1/r2)  ~  (r2vr> 


(II. D. 18) 

Returning  to  the  constitutive  law,  and  then  integrating  this  after  first 
substituting  eqs  17,18  yield 
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'dcff  t-u  f-vc) 


(II.D.19) 

Equations  13  and  19  define  the  rate  of  change  of  the  neck  radius  and 

grain  radius  respectively.  These  must  now  be  related  to  the  rate  of  change 

of  the  density  ratio.  If  V  is  the  volume  of  the  solid  ice  phase  in  a  cubical 

m 

region  containing  the  grain  of  radius  R  and  N  necks  of  length  L,  one  can  show 
that  and  the  volume  V  of  the  cubical  region  are 


m 


(2  /R2  -  A2  +  L)3 


(4-2N)7TR2/3  +  NttA2L/2 


+  (NTT/3)  (2R2  +  A2)  /R2-A2 


(II.D.20) 


(II.D.21) 


Vm  must  be  constant  if  the  matrix  material  is  incompressible.  Utili2ing  this 
constraint  gives  the  following  approximate  relationship  for  the  rate  of  change 
of  the  bond  length  L, 

£  ■  B2A  +  B2  R  (II. 0.22) 
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(II.D.23) 


where 

B,  =  -  ~r  {  (4  -  f  N)R2  +  N(2R2  +  A2)/ 3> 

2  N  A2  J 

Bx  -  “(2/NA2)  {nAL  +  2NRA/3  -  (2R2  +  A2) An/(3R) } (II.D.24) 

The  average  number  of  bonds  N  attached  to  a  grain  varies  with  density 
and  snow  type.  In  the  absence  of  any  substantial  data  on  the  variation  of 
N  with  density,  the  following  approximate  relation  is  assumed. 


N  -  4  (1  +  (p-0.30)/0.50)  (II.D.25) 

where  p  is  given  in  Mg  it  Very  little  data  are  available  to  verify  the 

accuracy  of  this  equation.  Kry  and  Gubler  made  studies  of  snow  structure, 

but  their  results  considered  a  relatively  narrow  density  range.  The  density 

ratio  can  also  be  expressed  as  the  ratio  V/V  .  Then,  the  rate  of  change  of 

m 

a  can  be  calculated  with  the  use  of  eqs  13,  14,  20,  21,  and  22.  The  result 
of  this  is 

a  =  A1  exp  (BjP)  +  A2  exp  (62p)  (II.D.26) 


where 


£  <r/a>2 

e2  =  |  [ (r/a) 2  f  -  i] 


(II .D. 27) 


Ax  =  (2KA/3A)  (-2A//  R2-  2  +  Bx)  exp  (-SQ/C) 


A2  =(KR/3A)  (2R//  R2-A2  +  B2)  exp  (-SQ/C) 


3(2 


/r2  -a2  +  l)2/v 


m 


Eq  26  describes  a  volumetric  behavior  for  snow  which  functionally  has  the 
correct  form.  However,  acceptable  accuracy  cannot  be  expected  since  factors 
such  as  work  hardening  of  the  matrix  material,  effects  of  random  grain  and 


neck  geometry,  and  intergranular  glide  have  not  been  accounted  for.  In  order 
to  adjust  this  equation  to  better  fit  experimental  data,  the  following  emperical 
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term  is  used 

C2  +  C3Po 

F  =  C.p^  (a/a)  (ll.D.28) 

l  o  o 

to  divide  the  pressure  p  in  eq  26. 

The  constitutive  equation  then  becomes 

a  -  exp  (61p/F)  +  A2  exp  (B2p/F)  (II.D.29) 

with  ,  A2 ,  6^  and  $2  remaining  unchanged  from  the  values  given  in  eqs. 

23,  24  and  27.  In  these  results,  all  coefficients  are  determined  by  the 
properties  of  ice  and  the  structure  of  snow,  except  for  the  three  constants 
C^,  C2  and  contained  in  eq.  28. 

11.  E  COMPARISON  WITH  EXPERIMENTAL  DATA  FOR  NECK  GROWTH  MODEL 

Abele  and  Gow  (1975,1976)  have  reported  results  of  their  studies  on 
the  high  rate  volumetric  properties  of  snow.  Their  1976  study  involved  the 
deformation  of  sifted  and  compacted  snow,  a  process  which  resulted  in 
initial  densities  higher  than  0.30  Mg  m  3.  Their  1975  work  was  concerned 
with  natural  undisturbed  snow  with  densities  ranging  from  0.1  Mg  m  3  to 
0.30  Mg  m  3  and  a  temperature  range  of  -3O°C<_0<  -2°C.  Their  experiments 
consisted  of  confined  compression  tests  and  have  already  been  described  in 
Section  II. C.  The  reader  is  referred  to  that  section  in  order  to  familiarize 
himself  with  the  test  procedures. 

Several  densities  were  considered  here  to  demonstrate  the  effect  of 
initial  density  on  the  stress  response  of  the  snow.  Generally,  the  average 
crystal  diameter  was  about  0.2  to  0.3  mm,  althouv  this  was  at  best  an 
estimate  obtained  from  the  thin  sections  illustrated  in  reports  by  Abele 
and  Gow.  This  snow  is  similar  to  fine  grained  seasonal  snow.  The  bond 
radius  was  taken  to  be  0.08  -  0.12  mm,  based  on  thin  sections  illustrated 
in  Abele 's  paper. 
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Figure  10  illustrates  a  comparison  of  theory  with  the  experimental 
data.  The  term  F  in  eqs  28  and  29  was  evaluated  by  determining  the  coefficients 
to  fit  eg  29  to  the  experimental  data  for  an  initial  density  of  0.30  Mg  m 
C^,  C2  and  were  found  to  be  0.112,  1.67,  and  6.40,  respectively.  As 
can  be  readily  seen,  for  all  densities  from  0.10  Mg  m  ^  to  0.60  Mg  m  3  the 
agreement  between  theory  and  experiment  is  excellent. 

Figure  11  illustrates  the  effect  of  bond  size  on  the  volumetric  properties 
of  snow.  In  Figure  11,  grain  size  is  held  constant  while  bond  size  is 
varied  over  a  range  characteristic  of  alpine  snow.  Figure  12  illustrates 
the  effect  of  deformation  rate  on  the  response  of  snow.  The  importance  of 
grain ssize  On. pressure  response  was  checked.  R  was  varied  over  an  order  of 
magnitude  while  the  ratio  A/R  was  kept  constant,  and  the  pressure  response 
was  found  to  vary  by  less  than  a  percentage  point,  so  that  one  may  conclude 
that  grain  size  is  a  second  order  effect. 

II. F  SIMPLIFIED  EQUATION  FOR  NECK  GROWTH  MODEL 

•  • 

Equations  13  and  19  can  be  used  to  calculate  A  and  R  during  volumetric 
compaction.  This  was  done  for  the  range  of  initial  densities  and  deformation 
rates  illustrated  in  Figures  10,  11  and  12,  and  R  was  found  to  be  at  least 
two  orders  of  magnitude  less  than  A.  This  is  due  to  the  fact  that  the 
grains  have  a  high  degree  of  lateral  constraint,  whereas  the  grain  bonds 
are  essentially  unrestrained.  This  result  allows  the  constitutive  equation 
to  be  simplified  by  neglecting  the  term  B2  exp  (SjFp) .  The  resulting  simplified 
constitutive  law  becomes 

.  B1P/F 

a  -  Axe  (Il.F.l) 
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and  this  equation  is  more  readily  inverted  to  obtain  p  as  a  function  of 
a  and  a.  This  equation  gives  results  virtually  identical  to  the  results 
illustrated  in  Figures  10  -  12. 

The  above  constitutive  equation  is  considered  to  be  valid  for  quazi- 
static  rates  of  loading.  When  very  large  rates,  such  as  those  associated 
with  stress  waves  are  incurred,  intergranular  dynamical  effects  such  as  those 
derived  in  the  pore  collapse  model  should  be  included.  This  results  with 
the  following  constitutive  law. 

2 


p  =  |^ln  (a/Ax)  +  V  ^  (a'a'S) 


^2.  (n.n/ri’i  (11. F. 2) 

In  summary,  a  rate  sensitive  volumetric  constitutive  law  based  on 
grain  and  grain  bond  deformation  has  been  developed  and  compared  with 
experimental  data.  As  can  be  seen  in  the  above,  the  developed  law  does 
accurately  describe  the  behavior  of  snow  for  a  given  deformation  rate  and 
for  a  range  of  initial  densities.  Aside  from  the  rheological  properties  of 
the  matrix  material,  one  of  the  most  important  parameters  is  the  ratio  A/r 
of  the  bond  and  grain  radii  in  the  undeformed  snow.  However,  the  importance 
of  this  parameter  becomes  less  significant  as  volumetric  strains  become 
large.  This  can  be  better  understood  by  first  considering  Figure  13,  which 
shows  the  variation  of  the  neck  radius  as  a  function  of  density  for  two 
different  intial  neck  radii.  It  can  be  seen  that  the  relative  differences 
in  bond  radii  decrease  during  the  compaction  process.  The  smaller  bonds 
experience  a  larger  stress  and  hence  flow  more  readily  than  the  large 
bonds,  thereby  growing  at  a  faster  rate. 

The  constitutive  law  was  found  to  lose  accuracy  at  the  higher  pressures 
where  the  theoretical  curves  are  terminated  in  Figure  10.  For  all  four 
cases  shown  in  Figure  10,  this  occurred  when  the  bond  radius  was  about  three 
fourths  the  value  of  the  grain  radius.  At  this  point,  th»»  adjacent  bonds 
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on  an  ice  grain  begin  to  interact  in  a  significant  degree  as  they  grow  into 
each  other.  The  kinematic  equations,  eq  22  -  27,  do  not  account  for  this, 
and  this  may  be  part  of  the  reason  for  the  loss  of  accuracy  at  the  higher 
densities. 

The  neck  growth  model  used  in  this  study  could  be  applied  to  other 
materials  such  as  porous  metals,  although  the  constitutive  law  for  the 
matrix  material  would  be  different  than  the  one  used  here.  However,  the 
kinematical  description  used  here  should  be  readily  adaptable  to  other 
porous  materials. 
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III  APPLICATION  TO  THE  PROBLEM  OF  VEHICLE  MOBILITY 


III. A  INTRODUCTORY  REMARKS 


The  problem  of  oversnow  mobility  can  generally  be  divided  into  two 
different  but  related  topics:  (1)  power  requirements  for  motion  in  snow 
and  (2)  traction  requirements.  The  latter  problem  has  probably  received 
■ore  attention,  since  sufficient  traction  and  drawbar  pull  capability  must 

be  developed  before  any  forward  motion  is  possible.  In  addition  the  problem 

i 

of  traction  capability  is  somewhat  simpler  to  define  mathematically,  at 
least  in  the  manner  that  this  problem  has  been  previously  treated  (Mellor 
1963) .  This  is  not  meant  to  imply  that  the  problem  of  traction  capability 
is  a  simple  problem,  since  the  stress  state  in  the  immediate  vicinity  of 
the  track  grousers  is  quite  complex.  Previous  studies  have  avoided  a 
detailed  stress  analysis  by  calculating  averaged  shear  stress  capability  of 
the  snow  at  the  track-snow  interface.  A  relation  similar  to  the  Coulomb- 
Mohr  criterion  is  used  to  find  the  shear  stress  in  terms  of  the  track 
pressure,  and  by  application  of  equilibrium,  the  tractive  capability  of  the 
vehicle  can  be  found.  This  approach  has  met  with  some  success  as  indicated 
by  Mellor  (1963) .  In  addition,  Harrison  (1956,  1975)  has  utilized  a  method 
used  for  soils  to  predict  drawbar  pull  capability,  but  he  has  pointed  out 
some  shortcomings  of  this  method. 

The  problem  of  calculation  of  vehicle  power  requirements  has  however 
not  been  studied  as  extensively  as  the  second  problem.  A  very  relevant 
reason  for  this  is  that  there  has  not  existed  a  theory  for  accurately, 
defining  the  volumetric  properties  of  snow  under  finite,  high  rate  com¬ 
paction  processes.  This  difficulty  has  in  part  been  relieved  by  the  development 
of  volumetric  constitutive  laws  such  as  presented  in  Section  II,  where  it 
was  shown  to  accurately  predict  the  response  of  medium-to-high  density  snow 


to  large  amplitude  high  rate  volumetric  deformations. 

In  the  case  of  a  shallow  snowpack,  the  pressure  bulb,  which  is  the 
region  of  compacted  snow  directly  beneath  the  vehicle  track,  extends  to  the 
ground.  This  enables  the  relatively  rigid  ground  to  lend  a  significant 
amount  of  support  to  the  vehicle.  If  the  snowpack  is  quite  shallow,  the 
effects  of  the  shear  stress  on  the  sides  of  the  pressure  bulb  may  be  neglected. 
The  pressure  inside  the  bulb  does  not  consequently  vary  signficantly  in  the 
vertical  direction,  thereby  allowing  the  rather  easy  calculation  of  power 
requirements  for  tracked  vehicles  in  snow. 

In  the  case  of  deep  snowpack,  the  pressure  bulb  receives  no  support 
from  the  ground,  and  the  bulb  must  be  supported  primarily  by  the  shear 
stresses  existing  on  the  bulb  wall.  As  a  consequence  the  pressure  bulb  can 
extend  a  significant  distance  into  the  snowpack  before  an  equilibrated 
state  is  reached.  Due  to  the  effects  of  the  shear  on  the  bulb  walls,  the 
pressure  and  the  density  distribution  in  the  wall  is  no  longer  uniform. 

The  constitutive  law  in  eq  II.B.47  is  used  to  estimate  energy  and 
power  requirements  for  over-snow  tracked  vehicles.  Since  high  rates  and 
large  amounts  of  compaction  are  involved,  constitutive  equations  developed 
by  other  researchers  are  generally  not  valid  for  this  case.  Mellor  (1974, 

1977)  gives  a  thorough  review  of  previous  constitutive  formulations  for 
snow,  and  one  can  readily  see  that  most  previous  work  is  valid  for  either 
small  strains  or  rates  well  below  those  that  can  be  generated  by  tracked 
vehicles  in  motion.  The  following  analysis,  therefore,  represents  one  of 
the  first  attempts  to  analytically  calcualte  power  requirements  for  tracked 
vehicles  which  generate  large  amounts  of  compaction  of  snow.  The  comparatively 
simple  problem  of  shallow  snowpack  is  first  studied,  and  then  deep  snowpack 
mobility  is  considered  in  the  following  section. 
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III.B  TRACKED  VEHICLE  MOBILITY  IN  SHALLOW  SNOWPACK  OF  MEDIUM  DENSITY 

Harrison*  (private  communication)  has  indicated  that  the  cross  section 
of  the  failure  zone  below  a  vehicle  track  takes  on  a  somewhat  rectangular 
shape.  Figure  14  gives  a  schematic  of  the  failure  region,  often  referred  to 
as  the  pressure  bulb.  It  is  assumed  that  shear  failure  occurs  along  the 
sides  of  the  failure  region,  and  that  the  material  directly  below  the 
vehicle  track  undergoes  unidirectional  compaction.  Outside  the  pressure 
bulb,  the  snow  undergoes  very  little  compaction  so  that  the  energy  dissipated 
would  be  negligible  when  compared  to  the  dissipation  occurring  within  the 
pressure  bulb.  The  depth  of  the  failure  region  depends  on  the  shear  strength 
of  the  snow,  and  some  of  the  energy  dissipation  is  invariably  due  to  deviatoric 
deformations.  However,  when  one  considers  the  massive  amount  of  compaction 
occurring  within  the  pressure  bulb,  neglecting  deviatoric  energy  dissipation 
becomes  a  valid  simplifying  approximation.  There  are  conditions  under 
which  such  an  assumption  may  not  be  reasonable.  Th*s  would  include  situations 
of  slippage  or  vehicle  turning.  But  if  a  vehicle  is  moving  along  a  straight 
path  and  is  not  experiencing  much  slipping,  the  above  approximation  should 
be  valid.  At  any  rate,  the  results  should  give  a  reasonable  lower  bound  on 
vehicle  power  requirements  and  allow  for  a  good  parametric  study  of  the 
problem.  The  following  analysis  would  be  particularly  valid  in  the  case  of 
shallow  snow  in  which  the  pressure  bulb  extends  to  the  ground. 

The  constitutive  equation  as  given  by  eq  II.B.47  can  be  used  to  calculate 
directly  the  energy  a  vehicle  expends  in  compacting  snow  during  its  passage. 
Stress  power  is  given  by 

W  =  ^  tr  (T  •  D)  (III.B. 1) 

a. 

* Research  Engineer,  U.S.  Army  Cold  Regions  Research  and  Engineering 
Laboratory,  Hanover,  New  Hampshire  03755. 
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where  T  and  D  are,  respectively,  the  Cauchy  stress  tensor  and  the  rate  of 
deformation  tensor,  tr(-)  is  the  trace  of  the  tensor  quantity  inside  the 
parenthesis.  The  stress  power,  W,  is  simply  the  rate,  per  unit  mass,  at 
which  work  is  being  done  internally  by  the  stresses,  ft  is  a  central  part 
of  the  first  and  second  laws  of  thermodynamics.  For  a  purely  viscous 
material  under  isothermal  conditions,  W  would  reduce  to  the  rate  of  energy 
disspiation.  For  a  viscoelastic  material,  the  stress  power  would  contribute 
to  both  the  rate  of  change  of  strain  energy  and  to  the  rate  of  energy 
dissipation. 

In  case  of  a  undirectional  deformation,  the  stress  power  associated 
with  compaction  is: 


ft  = 


£ 

p  3x 


(III.B.2) 


where  p  is  the  hydrostatic  pressure,  v  is  vertical  particle  velocity,  and  x 
is  the  vertical  coordinate  position  of  the  particle  during  deformation. 

The  continuity  equation  for  unidirectional  motion  is 


d£ 

dt 


+ 


0 


(III .B. 3) 


Eqs.  II.B.47,  2  and  3  can  then  be  used  to  obtain 


*  -  -  p  P  ° 
m 


(III.B.4) 


Integration  of  eq  4  results  with 

dt 


“-:T 


o  m 


(III  -B. 5) 

In  the  above  the  pressure  p  is  a  function  of  time.  Eq  II.B.47  may  be  in¬ 
verted  to  yield 

a  -  -  eQ(t'a)  (ill.B.6) 


where 


Q(t,a)  =  |  [F(t,a)  +  ln[a(a-l) ] ]  -  InA  (III.B.7) 

F(t,a)  =  |  {3ap(t)e(<f>a/aoV{Jln(^-))  -  2(SQ-C))  (III.B.8) 


The  pressure  loading  in  the  material  below  the  track  is  assumed  to 
have  .the  form 


p(t)  =  p*  —  [1  -  H(t-to)]  +  P*H(t-tQ) 


(1II.B.9) 


where  H(t)  is  the  Heavyside  step  function.  Figure  14  shows  the  nature  of  the 
pressure  distribution  under  the  track.  The  total  work  done  by  the  track  to 


a  unit  mass  of  snow  is  then 
*t* 

Jo  p 


w  £-P(t)  eQ(t'“dt 


(III. B. 10) 
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where  t*  is  the  duration  of  time  that  the  snow  is  under  the  track.  Eqs. 

7  to  10  can  then  be  used  to  study  the  energetics  of  oversnow  vehicle  travel. 

The  above  results  have  been  used  to  make  such  a  study.  By  substituting 
the  assumed  pressure  function  given  by  eq  9  into  eq  10  and  integrating  eq  10 
over  the  time  t*  corresponding  to  the  interval  that  an  element  of  snow  is 
under  the  vehicle  track,  the  actual  work  in  compressing  a  unit  volume  of 
undeformed  snow  is  found.  This  gives  a  direct  measure  of  track  efficiency 
for  a  given  set  of  parameters  such  as  vehicle  speed,  track  pressure,  track 
geometry,  and  snow  properties. 

Eq  2  gives  the  instantaneous  power/unit  mass  of  snow  while  the  material 
is  under  the  track,  and  this  expression  can  vary  considerably  during  the 
interval,  t*,  of  track  loading.  For  the  purpose  of  this  study,  the  average 
power  is  more  meaningful,  and  this  is  calculated  by  simply  dividing  the  total 
work  per  unit  mass,  W,  by  the  time  period,  t*,  required  to  produce  this  work. 


VEHICLE  SPEED  (m-s">) 


Therefore 


(III. B. 11) 


gives  the  power  per  unit  mass  is  subsequently  refrred  to  as  "specific  power". 
The  specific  power,  P,  would  have  to  be  augmented  by  a  factor  equal  to  the 
total  mass  of  snow  under  the  tracks,  if  one  wanted  to  find  the  total  power 
requirements. 

-Nominally,  a  track  length  of  5  meters  was  chosen.  Track  pressure  was 
4  -2 

limited  to  5x10  N-m  ,  which  admittedly  is  uncharacteristically  low  when 

considering  military  vehicles.  However,  such  pressures  are  more  realistic 

for  vehicles  such  as  snowmobiles.  Initial  snow  densities  studied  range 
-3  -3 

from  300  Ky-m  to  700  Kg-m  ,  and  vehicle  speeds  are  restricted  to  about 

ic  “1 

15  m-sec 

Figures  15  and  16  compare  specific  power  requirements  for  two  track 

loadings,  and  a  significant  difference  is  generated  when  one  increases  the 

4-2  4-2 

track  loading  from  1x10  N-m  “  to  5x10  N-m  .  In  particular,  the  efficiency 
of  the  lower  track  pressure  becomes  increasingly  apparent  at  higher  vehicle 
speeds.  This  is  due  to  decreasing  amount  of  snow  compaction  that  occurs 
while  the  snow  is  under  the  track  as  the  speed  increases.  At  the  lower 
track  loading  there  appears  to  exist  a  critical  speed  above  which  little 
increase  in  efficiency  is  achieved  with  higher  speeds.  For  the  higher 
pressure,  significant  amounts  of  compaction  continue  to  occur  at  the  higher 
vehicle  speeds.  These  high  pressures,  even  at  high  speeds,  forces  the 
vehicle  a  significant  distance  down  into  the  snowpack,  thereby  expending 
much  energy  in  snow  compaction .  At  lower  track  pressures,  the  compaction 
becomes  much  less  significant,  and  the  vehicle  tends  to  "ride-up"  on  the 
which  is  a  much  more  efficient  configuration.  The  relationship  between  the 
pressure  and  deformation  is  a  highly  nonlinear-  one,  as  evidenced  by  eq. 
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II. B. 47.  As  a  consequence,  one  can  expect  a  nonlinear  relationship  between 
specific  power  and  vehicle  speed. 

Notice  the  tremendous  increase  in  power  when  the  track  pressure  is 
4-2  4-2 

increased  from  1x10  N-m  to  5x10  N-m  .  Again,  this  can  be  attributed 
to  the  highly  nonlinear  relationship  between  p,  a,  and  a  in  the  constitutive 
equation . 

Figures  17  and  18  illustrate  the  variation  of  vehicle  power  with 

initial  snow  density,  tin  these  figures  the  relationship  between  specific 

power  and  density  is  illustrated  for  three  different  vehicle  velocities. 

As  can  readily  be  seen,  for  snow  with  initial  densities  above  300  Kg-m  3 ,  a 
4  -3 

track  loading  of  10  N-m  operates  fairly  efficiently.  However,  the  same 

4  -2 

cannot  be  said  for  a  track  loading  of  5x10  N-m  ,  where  good  efficiencies 
are  not  achieved  until  an  initial  snow  density  of  about  500  Kg-m  3  is 
reached . 

The  foregoing  calculations  indicate  that  much  can  be  done  in  an  analytical 
manner  to  make  parametric  studies  of  over-snow  mobility.  The  study  given 
here  considered  only  volumetric  effects,  and  as  such  would  represent  a 
lower  bound  on  total  energy  levels  absorbed  by  the  snow  when  compressed  by 
vehicle  tracks.  However,  this  estimation  should  be  good  if  slipping  is  not 
significant.  Central  to  such  an  analytical  study  is  the  availability  of  a 
constitutive  equation  which  can  accurately  represent  the  material  response 
to  large  compactions  at  high  rates. 

The  results  of  the  computations  illustrated  in  the  figures  show  an 
intricate  relationship  between  the  rate  at  which  energy  is  absorbed  by  the 
snow  and  such  parameters  as  track  pressure,  vehicle  speed,  and  initial  anow 
density.  Decreasing  the  track  pressure  substantially  reduces  the  snow 
density  required  to  allow  the  vehicle  to  move  efficiently  through  snow. 
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The  same  results  are  implied  for  the  speed  required  to  allow  the  vehicle  to 
begin  to  plane  out  onto  the  surface,  thereby  resulting  with  less  compaction 
of  snow.  The  figures  show  a  highly  nonlinear  relationship  between  these 
parameters  and  the  power  dissipation.  These  figures  also  indicate  there 
exists  critical  combinations  of  snow  density,  track  loading,  and  vehicle 
speed  which  allow  efficient  travel  over  snow. 

-The  calculations  were  made  by  assuming  a  uniform  compaction  of  the 

/ 

snow  directly  below  the  vehicle  track,  and  as  a  consequence  these  calculations 
would  be  particularly  accurate  for  shallow  snowpack,  in  which  pressure 
bulb  reaches  the  ground.  For  deep  snowpack,  compaction  within  the  pressure 
bulb  would  not  be  as  uniform,  and  the  results  would  have  to  be  interpreted 
in  terms  of  the  averaged  pressure  within  the  pressure  bulb.  However, 
without  knowing  the  shear  strength  of  the  snow  along  the  sides  of  the 
pressure  bulb,  these  calculations  would  be  difficult  to  arrive  at. 

Presently,  the  author  is  working  on  a  more  comprehensive  constitutive 
law  for  snow,  which,  when  completed,  should  provide  a  more  complete  description 
of  the  high  strain  rate  properties  of  snow.  In  addition  to  the  formulation 
of  a  more  comprehensive  constitutive  law,  much  more  experimental  data  is 
needed  so  that  the  variation  of  snow  properties  with  snow  type  can  be 
determined.  Density  alone  is  not  a  sufficient  parameter  for  doing  this, 
and  in-depth  studies  are  needed  to  determine  those  structural  parameters  of 
snow  which,  in  conjunction  with  density,  can  be  used  to  accurately  evaluate 
the  material  coefficients  to  make  calculations  such  as  done  here  more 
accurate.  When  a  set  of  constitutive  equations,  such  as  the  one  presented 
here,  are  available  for  a  wide  range  of  snow  densities  and  snow-types 
detailed  studies  can  then  be  made  on  vehicle  performance  in  snow  for  a  wide 
range  of  conditions.  Such  studies  can  form  guidelines  for  track  design. 
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III.C  TRACKED-VEHICLE  MOBILITY  IN  DEEP  SNOWPACK  OF  MEDIUM  DENSITY 

In  this  section,  we  assume  the  snowpack  is  so  deep  that  the  pressure 
bulb  does  not  extend  to  the  ground.  In  this  case,  the  pressure  bulb  must 
be  supported  by  the  surrounding  snowpack.  This  is  done  by  means  of  shear 
stresses  along  the  bulb  walls  and  a  normal  stress  along  the  bottom  of  the 
bulb. 

-Studies  by  Harrison  (1975)  have  indicated  that  in  deep  snowpack  the 
pressure  bulb  is  usually  very  close  to  a  rectangular  shape.  There  is  some 
spreading  of  the  bulb  below  the  track,  but  this  is  usually  not  very  significant. 
Figure  19  shows  schematically  a  typical  shape.  Often  local  inhomogeneous 
conditions  or  vehicle  turning  can  cause  the  vehicle  to  significantly  alter 
the  bulb  shape.  However,  if  the  vehicle  is  moving  in  a  straight  line,  the 
bulb  geometry  depicted  in  Figure  19a  is  realistic. 

Now  idealize  the  bulb  as  shown  in  Figure  19b.  The  hydrostatic  pressure, 
p  in  the  bulb  is  assumed  to  be  a  function  of  Y  and  t,  where  Y  is  the  verti¬ 
cal  Lagrangian  coordinate  of  a  particle  which  eventually  is  contained  in 
the  pressure  bulb  at  deformed  position  y  and  time  t.  The  shear  stresses 
along  the  wall  and  in  the  pressure  bulb  is  simply  the  critical  pressure  p 

c 

required  to  initiate  inelastic  compaction  of  the  snow.  Not  shown  at  the 
lower  surface  is  a  shear  stress  component,  which  must  be  zero  at  the  center 
of  the  bottom  surface  and  which  grows  as  the  two  lower  corners  are  approached. 
This  shear  stress  is  felt  to  partly  restrain  the  pressure  bulb  from  spreading, 
since  it  must  act  inward  toward  the  centerline  of  the  bulb. 

We  consider  now  the  Lagrangrian  forms  of  the  equilibrium  equation  and 
the  continuity  equation.  The  equilibrium  equation  and  continuity  equations 
are  respectively 
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Figure  19.  Actual  and  idealized  pressure  bulbs  below  vehicle  track. 
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and  T  are  respectively  the  normal  and  shear  components  of  the 

O 

Piola  stress  tensor  T.  The  Piola  stress  tensor  is  related  to  the  more 
familiar  Cauchy  stress  tensor,  T,  by  the  equation 


P°  -1 

T  =  —  F  .  T 

P 


(III.C.3) 


PQ  is  the  initial  density  of  snow,  and  F  is  the  inverse  of  the  deforma¬ 
tion  gradient  tensor.  Utilizing  eqs  3  and  1  and  taking  advantage  of  the 
uniformity  of  the  deformation  in  the  Z-direction  results  with  the  following 
reduced  form  of  the  equilibrium  equation: 


•V  2  p 

dp  __  O 

3y  p  w 
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(III.C.4) 


where  w  is  the  width  of  the  track,  and  T  is  shear  stress  on  the  bulb  wall. 


Integration  of  the  above  gives 


2  po  ry 

f,t)  =  -  — —  I  a 
p  w  J 


TdY '  +  p (0, t) 


(III.C.5) 


where  p(0,t)  is  the  pressure  produced  by  the  surface  loading. 
Eq  II.B.47  can  be  inverted  to  give 


rt  Q (t ' ,  p,  a) 

a  =  ~  J  e  dt'  +  a0 


(III.C.6) 


Q  =  j  IF(t'  a)  +  ln(o(a-l))]  -  In  A  (III.C.7) 

1 

P(t,  a)  =  ^PaepCy,  t)/(Jln  (^~ •)) -2  (So-C)  ]  (III.C.8) 


Equations  5  and  6  can  be  solved  simultaneously  to  yield  the  variation  of  a 
under  the  track.  In  order  to  solve  these  equations  the  shear  stress  T  must  be 
found.  In  his  earlier  study,  Mellor  (1963)  has  studied  the  shear  strength  of 
snow  and  have  shown  that  under  sufficiently  high  rates,  the  shear  stress  of  snow 
is  weakly  rate  dependent  but  depends  strongly  on  the  hydrostatic  pressure.  In 
view  of  these  results,  the  following  form  is  adopted  for  T. 

1  X  =  Tq  +  K  p  (III.C.9) 

Under  relatively  high  pressures,  the  portion  K  p  can  be  significantly  larger 
than  x  ,  so  that  Tq  can  be  neglected.  The  pressure  variation  then  becomes 

2  p  K  Ty  Tt  Q(a,  t') 

p (Y,  t)  =  +  -  /  /  e  dt*  p (Y 1 ,  t)dY'  +  p(0,t)  (III. C. 10) 

p  w  J  o  Jo 
m 

K  is  actually  a  function  of  the  velocity  of  the  particles  on  the  bulb  wall, 
but  here  it  will  be  assumed  constant,  since  the  rate  dependence  of  K  is 
not  well  defined. 

The  continuity  equation  can  be  put  in  the  form 


!1  0  0 


~  =  det  \l\  =  0  0 


M  3y  w 
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(III. C. 11) 


Integration  of  this  results  with 


(III. C. 12) 


y  is  simply  the  depth  of  the  rut  formed  by  the  tracks,  since  this  is  the 
o 

deformed  position  of  the  top  surface  Y  =  0.  At  the  bottom  of  the  pressure 

bulb,  y  =  Y,  since  there  is  no  significant  displacement  at  that  point. 

Therefore,  if  we  denote  Y  as  the  position  of  the  bottom  of  the  pressure  bulb, 
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eq  12  gives 


(1  -  — )  dY 
o 


(III. C. 13) 


For  a  given  problem,  eqs  6-13  can  be  solved  numerically  when  appropriate 

boundary  and  initial  conditions  are  specified.  Assuming  p  (0,t)  is 

known,  eq  II.B.48  is  first  evaluated  at  t=0  to  find  p  ,  the  critical 

c 

yield  pressure.  Eq  10  then  is  integrated  to  determine  the  value  of  Y  = 


Y_  at  which  p  =  p  .  Eq.  II.B.47  is  then  used  to  find  01  for  0<Y<Y„. 

Then  the  problem  is  step  forward  integrated  in  time  by  some  small  time 
increment  At,  and  the  procedure  is  repeated.  Doing  this  for  the  complete 
time  period  of  track  loading  gives  a  numerical  solution  defining  sinkage, 
energy  consumed,  a,  and  stress  distribution  and  bulb  depth. 

The  work  done  to  the  snow  can  easily  be  found  by  calculating  the  work 
expended  by  the  surface  pressure  in  compacting  the  snow.  At  the  upper 
surface,  the  work  is  given  by  the  expression 


dW  =  p(0, t)dy (0, t) 
The  work  rate  is  then 


(III. C. 14) 


dW 

dt 


p(0,t) 


SZ 

dt 


<0,t) 


(III. C. 15) 


This  is  the  work  rate/unit  time  per  unit  track  area.  Consequently  the  total 
vehicle  power  expended  in  compacting  the  snow  is 


-V  ■  £  f  g  ,0't’  at  B».e. 

J  O 

where  A„  is  the  area  of  the  vehicle  track  and  t*  is  the  time  that  the  snow 
T 

is  under  the  vehicle  track.  This  power  term  includes  the  work  dissipated 
due  to  shear  losses  along  the  bulb  wall  as  well  as  the  energy  dissipated  tj 
compaction  of  the  snow  in  the  bulb.  Not  reflected  would  be  energy  expended 
due  to  track  slippage  and  the  dissipation  of  energy  associated  with  the 
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large  deviatoric  deformations  in  the  immediate  vicinity  of  the  track  grousers. 

Idealize  now  the  track  pressure  by  the  distribution  shown  in  Figure 
20b.  The  pressure  is  assumed  to  build  up  linearly  to  a  constant  peak  pressure 
p*.  This  variation  is  given  by  the  equation 


P(0,t)  =  0 


t  <  t* 


p(0,t)  -  p  t  { x  —  H  (tt  ) ]  +  P*H (tt  ) ,  0  <  t  <  t* 
o  o  o  —  — 


(III. C. 17) 


H (t)  is  the  Heavyside  step  function,  and  tQ  is  the  entry  time,  i.e.  the  time 
increment  that  the  surface  snow  is  in  contact  with  the  front  part  of  the  track, 
t  is  determined  by  the  vehicle  sinkage  and  speed.  Figure  20b  illustrates  the 
temporal  variation  of  p(0,t).  Generally  the  track  pressure  has  a  periodic 
variation  determined  by  the  spacing  and  size  of  the  track  wheels.  For 
wheels  with  a  moderate  or  small  spacing,  this  periodic  fluctuation  may  be 
neglected. 

The  track  grousers  also  generate  stress  concentrations  within  the 
snowpack  in  a  localized  region  near  the  surface.  This  region  will  have  a 
very  complicated  stress  state  which  involves  both  large  normal  stresses  and 
large  deviatoric  stresses.  In  this  zone,  a  significant  amount  of  energy  is 
dissipated  by  both  shearing  and  compaction,  and  the  exact  nature  of  the 
deformation  within  this  area  would  be  determined  by  grouser  size,  geometry 
and  spacing  as  well  as  the  amount  of  slippage  and  nominal  track  pressure. 

The  analysis  made  in  this  section  does  not  have  the  capability  of 
calculating  the  energy  absorbed  in  this  zone  within  the  localized  effect  of 
the  track  grousers.  Here  we  find  the  energy  absorbed  through  compaction  of 
the  snow  in  the  pressure  bulb  and  the  shearing  deformation  within  the 
pressure  bulb  and  along  the  bulb  walls.  As  a  consequence  these  calculations 
provide  a  lower  bound  on  the  total  energy  absorbed  by  the  snow.  In  the 
absence  of  vehicle  slipping,  the  shear  zone  in  the  grouser  region  may  be 


neglected  in  comparison  to  the  much  larger  pressure  bulb. 
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Figure  22.  Effect  of  initial  snow  density  on  power  consumption 
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Vehicle  pi'rfdriMiio;  ran  1...  d**t.ermi  tied  ;  ■  it aim -t.r  j  ca  1  1  y  in  terms  of  a 
number  of  pur.un*.  tor:;,  i  nclud  i  ng  track  geometry,  tr  ack  prossur- • ,  vehi  c ! 
s:  ood,  and  six >w  densi  ty.  By  track  geometry,  we  mean  track  width,  length, 

•  nil",  b  agt !i  and  augl-  of  attack. 

In  adui  ti-.m  t  a  calculatin';  vehicle  j")/.'  i‘  requirements,  :;uci.  i  n format  ion 
as  it  i  nkuqe  ,  kill  b  de\  tii,  and  density  and  ;  tossur**  ;  rotiln  within  tin.-  r  u ! 
can  is-  found,  Not  all.  of  this  will  la*.-  dom  here ,  since  this  'would  require 
to*.  much  stile.:,  and  the  main  interest,  is  the  evaluation  of  power  requirements. 

The  S  usie  track  geometry  and  vehicle  we  iqht  chosen  for  the  calculations 
a •  •  those  of  the  ordinance  M3A4  high-speed  tractor.  This  vehicle  was 
ori'iinally  designed  for  transporting  personnel  and  light  cargo  over  soft 
t«  rr.iin.  The  Ml  A  4  has  a  33f>  horsepower  gasoline  engine,  weighs  about 
t Ive  tons,  and  has  a  track  length  of  afoul  3.0  meters  and  a  track  width 
>*’  .i.Jl  me  tors .  The  nominal  track  pressure  varies  from  0.43  bar  upward, 
depending  on  the  load .  Consequently  the  M5A4  does  not  have  good  capability 
mu  ej  •  ■;  ,i  t  i  no  in  de.-p  seasonal  snowpuek  but  does  have  a  reasonable  mobility 
i  r:  per*  nnial  r.nowpack  such  as  on  the  Oroenland  icecap  or  on  medium  density 
■  .■  ■  a  son.  1 1  snow’,  generally  those  densities  ar*-  in  excess  of  3iiy  K.g-m  J . 

1  i'lui'e  .d  -37  summai  i  xo  the  results  of  this  study.  The  results  are 
iii  own  in  tel  ms  of  horsepower,  since  this  is  the  most,  recognizable  unit  cf 
j  (<wer.  The  track  pressures  considered  were  for  the  most  part  larger  than 
bars,  since  for  medium-to-high  d*  -nsi  ty  snow,  track  pressures  below  0..! 
bat :  result  with  insignificant  energy  levels  consumed  in  the  pressure  bulb. 

iiguies  31.  and  33  illustrate  the  very  strong  dependence*  of  energy 
consumption  on  initial  density  and  nominal  track  pressure.  For  instance  by 
tm-iely  increasing  the  pressure  from  0.3  liars  to  0.9  bars,  an  increase  of  an 
order  of  magnitude  in  energy  consumptin  results  for  snow  with  a  density  of 
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TRACK  LENGTH = 3.0m 
TRACK  WIDTH = 0.31  m 


300  Kg-m  \  Conversely  snow  with  an  initial  density  of  500  Kg-m  consumes 

only  about  10%  as  much  energy  as  snow  with  an  initial  density  of  300  Kg-m 

as  shown  in  Figure  5.  As  indicated  by  Brown  (1978b)  and  evidenced  in 

Figure  22,  for  a  given  pressure  P*,  there  exists  a  critical  density  above 

which  very  little  compaction  and  energy  consumption  takes  place.  Figure  22 

_  3 

shows  that  for  a  density  above  450  Kg-m  a  track  pressure  of  0.6  bar 
produces  very  little  compaction. 

Figure  23  demonstrates  the  effect  of  track  geometry  on  track  efficiency. 
The  M5A4  has  a  length-width  ratio  of  about  10,  but  Figure  23  shows  how 
energy  consumption  would  vary  for  a  L/w  variation  from  2  up  to  10.  This 
figure  shows  the  obvious  advantage  of  a  long  narrow  track.  The  reason  for 
this  dependence  lies  in  part  with  the  role  that  the  shear  stress  plays  in 
supporting  the  pressure  bulb.  One  can  readily  see  by  Figure  19  that  for  a 
wide  track,  the  pressure  bulb  would  have  to  be  deeper  in  order  for  the  wall 
shear  stresses  to  support  the  bulb.  This  can  also  be  seen  by  noting  eq  4 
in  which  the  track  width  w  inversely  affects  the  pressure  gradient  in  the 
bulb.  An  increased  w  decreases  9p/3y,  thereby  increasing  the  depth  to 
which  the  bulb  penetrates  before  the  bulb  pressure  decreases  to  the  critical 
yield  pressure  pc* 

Figure  24  shows  the  variation  of  bulb  depth  with  track  pressure  for 
three  initial  snow  densities.  As  should  be  expected,  bulb  depth  decreases 
with  decreasing  pressure  and  increasing  density.  The  variation  of  bulb 
depth  with  track  length-width  ratio  is  illustrated  in  Figure  25,  and  Figure 
26  gives  the  vertical  variation  of  bulb  pressure  for  one  particular  case. 

As  can  be  seen,  the  dpeth  is  critically  dependent  on  the  value  of  pc,  the 
yield  pressure  given  by  eq.  II.B.48.  Over  the  bottom  0.4  meters  of  the 
bulb,  the  pressure  p  goes  through  a  relatively  small  variation,  so  changing 


73 


MILES  PER  HOUR 


Figure  27.  Variation  of  compactive  force  with  vehicle  speed 
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by  50  to  100%  could  case  a  significant  change  in  the  bulb  depth.  This, 
however  would  not  affect  power  consumption  much,  since  most  of  the  energy 
is  consumed  in  the  upper  portion  of  the  bulb. 

Finally,  Figure  27  shows  the  compactive  force  generated  by  the  motion 
of  the  vehicle  in  the  snow.  This  was  obtained  by  integrating  the  horizontal 
component  of  the  track  pressure  over  the  track  area.  Note  here  the  improvement 
of  track  efficiency  with  vehicle  speed. 

The  results  of  this  study  show  that  the  method  described  here  has  the 
potential  of  predicting  power  requirements  of  oversnow  vehicles  for  a  wide 
range  of  conditions.  For  example  the  formulation  developed  here  enables 
one  to  make  a  parametric  study  of  the  effect  of  a  large  number  of  factors 
on  energy  consumption  rate.  This  includes  such  factors  as  nominal  track 
pressure,  length-width  ratio,  snow  properties,  angle  of  attack,  etc. 

Therefore  this  formulation  could  be  of  some  use  as  an  analysis  and  design 
tool. 

The  calculations  used  in  the  example  were  for  snow  with  a  uniform 
density  and  a  track  which  develops  a  uniform  track  pressure.  Such  an 
assumption  is  not  necessary,  and  a  stratified  snowpack  and  nonconstant 
track  pressure  could  have  been  used.  However,  to  demonstrate  the  formulation, 
the  simpler  problem  was  considered  here. 

Presently  there  is  not  much  data  available  for  comparison  with  the 
example  given  here.  Harrison  (1975)  has  detailed  results  for  several 
vehicles  in  seasonal  snowpack  in  which  th  snow  density  in  the  upper  50  cm 
were  less  than  300  Kg-m  which  makes  direct  comparison  with  the  example 
given  here  difficult.  The  constitutive  law  given  in  eq  II.B.47  is  not 

-3 

considered  to  be  accurate  for  densities  less  than  300  Kg-m  ,  so  cons¬ 
equently  an  example  for  lower  density  snow  was  not  possible. 
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The  general  behavior  of  the  results  demonstrated  in  Figures  21-27 
agree  with  what  has  been  observed  in  field  studies  reported  by  Harrison 
(1975) .  Whether  or  not  there  is  a  good  quantitative  agreement  will  be 
determined  when  either  new  field  data  is  available  or  when  the  model  developed 
here  has  been  generalized  to  include  conditions  of  lower  density  snow. 

What  has  been  demonstrated,  however,  is  the  apparent  usefulness  of  this 

model  for  design  and  analysis  purposes. 

I 

1II.D  VEHICLE  PERFORMANCE  UNDER  DIFFERENT  SNOWPACK  CONDITIONS. 

In  the  previous  two  sections,  vehicle  performance  was  evaluated  for 
dry  snowpack  with  densities  greater  than  300  Kg  m  3.  The  pore  collapse 
equation  was  used  for  these  calculations,  since  it  is  valid  for  this  type 
of  snow  and  is  relatively  easy  to  use. 

In  this  section  the  neck  growth  model  is  used.  This  allows  an  evalua¬ 
tion  of  vehicle  performance  in  low  density  snow  as  well  as  higher  density 
snow.  In  addition,  since  grain  size  and  bond  diameter  can  be  varied,  the 
effects  of  snow  type  can  be  studied. 

Normally  a  vehicle  will  have  to  be  operated  under  a  variety  of  snow¬ 
pack  conditions.  For  instance  during  mid  winter,  operation  on  low  density 
dry  snowpack  may  be  required,  whereas  during  the  spring,  unsaturated  wet 
snow  will  be  incurred.  It  would  be  useful  to  be  able  to  predict  how  much  • 
variations  in  snowpack  conditions  affect  vehicle  performance. 

Recently  Abele  (1981)  has  investigated  the  volumetric  properties  of 
wet  snow.  Of  particular  interest  is  his  work  on  unsaturated  snow  (free 
water  content  less  than  8%  by  volume) .  A  volumetric  constitutive  law  for 
this  snow  (Brown,  unpublished)  has  been  developed  and  could  be  used  in  a 
number  of  problems  of  interest.  This  constitutive  equation  has  a  form 
similar  to  the  poro  collapse  model  developed  earlier,  i.e. 


i 

V 

1 


ft 


p  „  g.~<fa/ao  [2<so  -  C)  4  C  m  ,]  (III. D.  1) 

which  is  essentially  the  same  equation  as  was  used  to  describe  the  behavior  of 
high  density  snow  (Brown,  1979) .  Abele  (1980)  indicates  that  wet  snow  is  somewhat 
more  rate  dependent  than  dry  snow,  and  that  the  material  is  somewhat  softer 
than  dry  snow  at  the  same  density.  Also  the  work  hardening  term  <j>  was 
found  to  be  dependent  on  the  free  water  content.  Utilizing  Abele' s  experimental 
results,  eq  11  was  used  to  approximate  the  properties  of  wet  snow  with 
essentially  no  free  water  and  wet  snow  with  7%  free  water  by  weight.  It 
should  be  noted  that  to  date,  Abele1 s  data  is  the  only  high  rate  test  data 
for  wet  snow,  so  the  constitutive  equation  given  by  eq  11  cannot  be  expected 
to  be  extremely  accurate.  There  just  is  not  yet  enough  data  to  thoroughly 
verify  the  accuracy  of  any  equation.  However,  the  results  found  here  for 
wet  snow  should  give  a  good  qualitative  evaluation  of  mobility  in  wet  snow. 

A  procedure  essentially  the  same  as  used  in  the  earlier  two  sections 
can  be  used  to  analyze  vehicle  performance,  so  only  the  results  are  given 
here.  The  same  track  loading  is  assumed  here  as  was  assumed  in  Section 
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III.C. 

For  dry  snow,  a  range  of  densities  from  200  Kg  m  3  to  500  Kg  m  3  were 
evaluated  for  speeds  ranging  from  1  ms  1  to  10  ms  1.  Also  the  effect  of 
track  pressure  on  vehicle  power  requirements  was  also  investigated.  These 
results  are  all  illustrated  in  Figures  28-32.  The  unit  of  power  used  here 
is  horsepower  as  it  is  more  familiar  to  people.  Figure  31  illustrates  the 
effect  of  snow  type  on  vehicle  performance  in  snow.  Plotted  in  that  figure 
is  the  variation  of  vehicle  horsepower  with  the  ratio  A/R.  For  well  sintered 
mid-season  snowpack,  the  mean  radius,  r,  is  on  the  order  of  0.1  mm,  whereas 
A  may  average  as  large  as  0.04  rnir.  However,  snow  which  has  been  under 
'•-cmperature  gradient  effects  will  metamorphose  until  the  grain  size  R  may 
be  as  large  as  2  mm. 
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A  relative  comparison  of  vehicle  power  requirements  in  wet  snow  and 

dry  snow  is  illustrated  in  Figure  32  where  results  are  shown  for  wet  snow 

with  0%  free  water  and  7%  free  water.  A  density  of  400  Kg  m  ^  was  chosen 

-3 

since  late  season  wet  snow  usually  has  a  density  in  excess  of  400  Kg  m 
For  the  wet  snow,  eq  1  was  used  in  evaluating  the  vehicle  performance. 

Figure  28  shows  the  dramatic  increase  in  energy  requirements  with 
decreasing  snowpack  density.  The  sinkage  increases  significantly,  and  this 
additional  sinkage  increases  the  thrust  required  to  push  the  snow  just 
ahead  of  the  track  downward  and  forward  as  the  track  passes  over  the  snowpack. 
Comparisons  in  Figure  28  are  for  two  vehicle  speeds  while  such  parameters 
as  intergranular  bonding  were  kept  constant. 

One  additional  point  should  be  made  here.  As  rut  depth  increases,  the 
forward  thrust  applied  by  the  tracks  also  must  increase.  This  results  with 
an  increased  shearing  force  applied  to  the  snowpack  in  a  direction  parallel 
to  the  direction  of  vehicle  travel.  This  increased  shearing  force  results 
with  additional  horsepower  consumption  within  the  snowpack  and  would  have 
to  be  added  to  the  numbers  arrived  at  here  since  shearing  effects  were  not 
included  in  the  analysis.  Consequently,  the  solution  found  here  should 
represent  a  lower  bound  on  the  actual  energy  consumed  within  the  snowpack. 
Studies  with  w.  Harrison  of  the  U.  S.  Army  Cold  Regions  Research  and 
Engineering  Laboratory  show  that  shearing  effects  are  not  significant  as 
long  as  the  vehicle  is  not  experiencing  slippage. 

Figures  29  and  30  show  that  power  consumption  increases  rapidly  with 
track  pressure  but  that  when  a  critical  vehicle  speed  is  reached,  vehicle 
travel  becomes  relatively  efficient.  This  critical  speed  depends  on  snowpack 
density,  track  pressure  and  intergranular  bonding. 


Finally  the  variation  of  vehicle  power  required  for  travel  in  temperature 


gradient  snow  increases  dramatically  as  the  temperature  gradient  effects 
become  more  significant.  Mid  season  equi-temperature  snow  is  the  pre¬ 
dominant  snow  type  for  normal  snowpack.  For  instance  in  a  two  meter  deep 
snowpack,  a  full  .5  meters  may  consist  of  this  snow  type.  Grain  size  is 
normally  0.1  -  0.3  mm  in  diameter,  and  intergranular  bonding  is  well  estab¬ 
lished,  the  ratio  A/R  often  being  on  the  order  of  0.3  to  0.5.  However, 
under  the  influence  of  a  temperature  gradient,  grain  size  increases.  This 
rate  of  increase  depends  on  temperature,  so  that  the  metamorphism  proceeds 
most  quickly  near  the  bottom  of  the  snowpack  where  the  temperatures  are 
normally  above  -5°C.  Quite  often  in  a  one  meter  early  season  snowpack 
which  experiences  a  cold  period,  the  entire  snowpack  can  show  these  tempera¬ 
ture  gradient  effects.  It  is  certainly  possible  that  the  ratio  A/r  could 
decrease  to  values  less  than  0.1  to  produce  very  weak,  collapsing  snowpack. 
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Figure  32.  Comparison  of  power  requirements  for  wet  snow  and  dry  snow 
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which  is  density  dependent.  The  pressure  term  is  a  function  of  a  and 

« 

a  and  is  given  by  eq  II.B.47.  For  the  settlement  problem,  this  formulation 

-4  -1 

was  found  to  work  well.  Generally  strain  rates  less  than  10  s  were 
incurred  in  the  problem.  Whether  or  not  it  would  work  at  all  well  for 


higher  strain  rates  associated  with  vehicle  mobility  has  not  been  determined. 
Also  the  constitutive  equation  is  somewhat  complicated,  and  it  should  be 
possible  to  formulate  a  simpler  constitutive  equation  for  the  vehicle 
mobility  problem.  At  any  rate,  if  an  accurate  representation  of  the  deformation 
in  the  bulb  below  the  track  is  to  be  determined  for  conditions  where  vehicle 
slippage  is  occurring,  a  more  general  equation  than  used  in  this  paper  is 
needed.  Assuming  an  acceptable  constitutive  relation  for  snow  can  be 
found,  solution  procedures  similar  to  that  of  Yong  et  al  (1978)  can  be 
performed.  Yong's  work  to  date  represents  the  best  numerical  work  on  this 
type  of  problem.  Their  approach  would  have  to  be  modified  to  account  for 
large  snowpack  strains  and  somewhat  more  complicated  material  properties, 
but  their  finite  element  approach  should  give  useable  results. 
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XV  APPLICATION  TO  SHOCKWAVE  PROPAGATION 
IV. A  INTRODUCTORY  REMARKS 

Stress  waves  in  snow  are  a  problem  that  has  not  been  studied  as  extensively 
as  other  subjects.  Practically  all  of  the  previous  work  on  the  mechanical 
properties  of  snow  has  been  restricted  to  quasi-static  conditions  in  which 
inertial  effects  are  negligible.  Napadensky  (1964)  first  investigated  the 
dynamic  properties  of  snow  under  the  effects  of  shock  waves.  More  recently, 
Johnson  (1978)  considered  in  detail  the  propagation  of  elastic  sonic  waves 
in  snow.  Wakahama  and  Sato  (1977)  Wisotski  and  Snyder  (1964) ,  and  Gubler 
(1977)  all  conducted  experimental  investigations  of  stress  waves  in  snow. 

Mellor  (1977)  reviewed  previous  work  on  shock  waves  in  snow,  as  did  Johnson 
(1978) .  But  to  date  virtually  no  work  performed  has  utilized  a  material 
constitutive  equation  to  investigate  the  propagation  of  inelastic  stress 
waves  in  snow.  This  is  not  surprising,  since  a  valid  constitutive  law  for 
snow  has  been  lacking,  and  consequently  previous  studies  have  been  restricted 
to  the  use  of  mass  and  momentum  balance  principles.  However,  o..ce  a  constitutive 
law  is  found  to  accurately  describe  the  behavior  of  snow  under  rates  of 
loading  characteristics  of  stress  waves,  a  much  more  detailed  analysis  of 
stress  waves  can  be  made. 

A  number  of  problems  require  a  detailed  knowledge  of  the  response  of 
snow  to  shock  waves.  The  relative  effectiveness  of  in-snow  and  airborne 
explosives  for  initiating  avalanches  is  one  such  example.  To  date  the 
relative  effectiveness  of  explosives  detonated  in  the  air,  or  in  snovpack, 
or  on  the  ground  has  not  been  determined.  Avalanche  experts  still  do  not 
agree  on  optimum  explosive  speed  or  charge  size.  Another  case  in  which 
stress  waves  are  generated  is  projectile  impact.  A  related  problem  is 
avalanche  impact  on  structures.  In  all  of  these  problems,  solution  of  the 
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usual  balance  principles  can  yield  some  information;  but  without  an  appropriate 
constitutive  law,  such  questions  concerning  stress  wave  attenuation  or 
alteration  of  stress  wave  profile  cannot  be  answered. 

At  this  point  some  definitions  are  in  order.  An  elastic “wave  is  one 
in  which  no  internal  dissipation  takes  place.  The  material  responds  elastically 
to  the  stress  wave,  and  no  attenuation  can  be  attributed  to  inelastic 
deformation  of  the  material.  A  plastic  wave  is  one  in  which  material 
yielding  and  viscous  flow  occur.  This  is  often  referred  to  as  a  nonlinear 
wave.  This  is  an  inexact  definition,  however,  since  material  nonlinearity 
cam  occur  in  the  absence  of  material  yielding  and  plastic  flow.  However, 
since  at  high  deformation  rates  snow  remains  practically  linear  up  to  the 
point  where  yielding  begins,  any  differentiation  between  plastic  waves  and 
nonlinear  waves  is  meaningless. 

In  stress  waves,  the  jump  in  a  variable  (such  as  strain,  pressure,  and 
energy)  is  defined  as  the  difference  between  the  values  of  the  variable 
just  in  front  of  and  behind  the  wave  front.  Much  of  the  previous  analytical 
work  on  stress  waves  has  been  restricted  to  the  determination  of  these 
jumps.  A  shock  wave  is  a  wave  in  which  the  displacement  of  a  particle  is 
continuous  across  a  wave  front  but  in  which  the  particle  velocity  experiences  i 

a  jump.  An  acceleration  wave  is  one  in  which  the  particle  velocity  and  f 

displacement  are  continuous  across  a  wave  front,  but  in  which  the  accelera¬ 
tion  has  a  jump  across  the  wave.  The  surface  representing  a  wave  front  is 
often  referred  to  as  a  singular  surface,  since  discontinuities  in  acceleration 
and/or  velocity  can  occur  at  this  surface. 

A  steady  wave  is  defined  here  as  one  in  which  the  wave  speed  propagates 
at  a  constant  speed  V.  One  can  show  that  the  wave  amplitude  remains  constant 
for  such  a  wave.  Even  in  plastic  waves  such  a  condition  can  be  established, 
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although  generally  this  condition  is  short-lived.  Steady  waves  have  been 


studied  extensively  since  they  are  mathematically  easy  to  investigate. 

Nonsteady  waves  are  simply  those  waves  that  lack  the  characteristics  of 
steady  waves;  i.e.,  their  amplitudes  may  either  grow  or  attenuate. 

IV. B  GOVERNING  EQUATIONS  FOR  PLANE  SHOCKWAVES 

In  the  following  sections,  differentiation  of  the  constitutive  law  will  be 


necessary.  This  can  be  done  in  terms  of  the  moduli: 

_-3£  (a, a, ci) 

fcT  ~  3a 


E1  =  la 


E  = 

2  3a 


(IV. B. 1) 


(IV. B. 2) 


(IV. B. 3) 


Et,  and  E2  are,  respectively,  the  tangent,  rate,  and  acceleration  moduli. 

The  derivatives  3p/3x  and  dp/dt  then  can  be  shown  to  have  the  forms : 


3p  3a  ,  _  3d  „  3a 

3X  ='ET  W.  +  E1  H  +  E2  5x 


to  =  -ETa+Ela+;:2a 


(IV. B. 4) 


(IV. B. 5) 


where  X  is  a  coordinate  variable. 

Consider  now  balance  laws  for  stress  waves.  We  use  here  strictly  a 
mechanical  theory  and  do  not  study  restrictions  invoked  by  the  energy  law 
or  the  second  law  of  thermodynamics.  The  two  laws  we  consider  are  those  of 
mass  balance  and  momentum  balance. 

The  momentum  balance  for  a  plane  wave  propagating  in  the  X  direction 
requires 
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where  is  the  body  force  acting  in  the  X  direction  and  V  is  the  particle  velocity. 

If  the  state  of  stress  is  dominated  by  the  pressure  p,  and  if  body  forces  are 
negligible,  eq  6  can  be  reduced  to: 


3v 


3x  K0  3t 

The  mass  balance  equation  (or  continuity  equation)  is 


(IV. B. 7) 


p0  3p  3v 
p  3t  +  3x 
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Since  a  =  p  /  p  =  p_a  /p,  we  can  get 
m  u  u 


3v  1  3a 
3X  _  aQ  3t 


(IV. B. 8) 


(IV. B. 9) 


as  the  form  of  the  mass  balance  equation  which  is  used  later. 

Consider  a  one-dimensional  stress  wave  propagating  through  a  medium 
such  as  snow.  We  define  as  a  wave  (or  wave  front)  a  smooth  one-parameter 
family  of  points  Y  (t)  ,  °°,<  t<  such  that  Y  (t)  gives  the  material  point 
(or  particle  X)  at  which  the  wave  is  located  at  time  t.  X  is  the  position 
of  a  particle  in  the  reference  configuration,  which  here  will  be  the  un¬ 
deformed  configuration,  x  =  x(X,t)  is  the  position  of  a  particle  X  at  time 
t  and  is  therefore  the  deformed  position.  The  intrinsic  velocity  of  the 
wave  is 


V  =  ft  Y(t)  (IV. B . 10) 

which  is  the  velocity  of  propagation  relative  to  the  undeformed  position  of 
the  body. 
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I  Let  f  be  any  variable,  say  density  or  stress.  Assume  f(X,t)  is  a 

|  function  of  position  X  and  time  t  and  is  of  class  in  X  and  t  except  at  X 

=  Y,  where  f  has  a  jump  on  discontinuity;  i.e.,  the  values  of  f  just  in 
|  front  of  and  behind  the  wave  front  (X  =  Y)  have  different  values.  We 

denote  this  jump  by  the  expression  [f] ,  or 

!  [f]  =  f-f+  (IV.B.ll) 

■  i 

|  where  j 

f  =  lim  f  (X,  t) 

|  X  -*•  Y  (IV.B.12) 

X  <  Y 

f+  =  lim  f (X,t) 

X  Y  (IV .B.  13) 

X  >  Y. 

+  — 

f  and  f  represent,  respectively,  the  values  of  f  just  in  front  of  and 
just  behind  the  wave  front.  For  instance,  in  a  shock  wave  we  would  have 
(v)y  0,  [v]  ^  0,  and  [x]  =  0;  i.e.,  the  particle  acceleration  would  undergo 
a  jump  across  the  shock  wave,  but  particle  position  would  not.  The  above 
equations  were  all  defined  with  the  assumption  that  V  >  0,  so  that  the 
wave  is  propagating  in  the  positive  coordinate  direction. 

Now,  consider  derivatives  of  [fj,  since  these  will  be  used  later  in 
the  analysis  of  shock  waves.  Assume  f(X,t)  has  a  jump  discontinuity  at  the 
wave  front,  X  =  Y (t) ,  but  otherwise  is  continuous  and  differentiable. 

Thus,  clearly  the  jump  [f]  is  a  function  of  time  only  through  the  position 
X  =  Y (t) .  The  derivative  of  I f J  follows  from  the  definition  of  the  jump: 

dt  |f)  “  at  (  ,iY 


(t),  t)-f(Y  (t) 


„) 


92 


I 

I 


3f(Y  ,t)  dY  3f 

ax  at  at 


df(Y  ,t)  dY 

ax  dt 


3f 

at 


from  which  we  get 


df  lfl 
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»  m 

(IV.B.14) 


However,  if  we  assume  f  is  continuous  across  the  wave  front,  then  If]  =  0, 
and  eq  14  gives 


v  [14  -  -  [it] 


(IV.B.15) 


Equation  15  places  restrictions  on  the  jumps  in  the  derivatives  of  f  if  f 
is  continuous  at  Y.  Equation  15  is  called  Maxwell's  theorem,  although  both 
eq  14  and  15  are  also  often  referred  to  collectively  as  the  compatibility 
conditions  for  singular  surfaces. 

The  compatibility  equations  can  now  be  applied  to  the  balance  laws. 
Recalling  that  a  shock  wave  generates  jump  discontinuities  in  v,v,  and  3v/3x, 
whereas  x  isccontinuous  across  the  wave.  Maxwell's  theorem  gives 


IV]  =-V  (IV.B.16) 

To  calculate  the  pressure  jump  across  a  pressure  wave,  integrate  eq  7  from 

X^  to  X„,  where  X„  <  Y(t)  <  XD,  to  obtain 
a  p  a  p 


-  p(Xg,t)tp(Xa,t)  =  fx  B  P0  X 


x  dX 

Or,  by  breaking  the  integral  on  the  rightnhand;;side\into;.twopparts,  we  get 

-p(Xg,t)  +  p<Xa,  t)  = 

'y(t) 

K 

a 


fe  (rv-x;6,  oo  h 


(IV.B.17) 
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Now,  taking  the  limits  of  Y~(t),  X„  Y  (t) ,  and  using  eq  16  and 

a  p 

Leibnitz's  rule  to  differentiate  the  integral  on  the  right-hand  side,  we  obtain 


Ip]  ■  P0  V  lx]  =  pQ  V  lv] 


(IV.B.18) 


This  equation  is  a  familiar  jump  equation  that  relates  the  pressure  jump 
across  the  wave  to  the  jump  in  the  particle  velocity  and  wave  speed.  A 
similar  procedure  applied  to  eq  9  results  in 


(.1  t  ^  Id 


(IV.B.19) 


Therefore,  we  see  that  the  material  has  a  jump  in  a  if  the  wave  is  a  shock 
wave,  since  [v]  is  nonzero.  Combining  eq  18  and  19  gives 


(1V.B.20) 


which  is  a  well  known  relationship  for  shock  waves  relating  pressure  jump 
across  a  wave  to  the  jump  in  density  ratio. 

Evaluating  eq  7  both  ahead  Of  and  behind  the  wave  front,  then  letting 
X  approach  Y(t)  and  subtracting  these  two  forms  of  eq  7  lead  to: 


[ft]  =  “PQ  ™ 


(IV.B.21) 


Equations  18-21  are  all  widely  recognized  jump  equations  for  stress  waves 
and  are  applicable  to  both  steady  and  nonsteady  waves. 

We  now  analyze  in  detail  the  propagation  of  steady  waves  before  considering 
the  details  of  nonsteady  stress  waves. 

IV. C.  STEADY  SHOCK  WAVES  IN  SNOW 

To  investigate  further  the  properties  of  steady  waves  we  now  return  to 
the  jump  equations,  eq  IV. B. 18-21.  If  the  material  is  at  rest  just  prior 
to  the  stress  wave  arrival,  the  jumps  [p]  and  [v]  equal  the  values  of  p  and 
v  as  the  wave  passes.  Denoting  the  steady  values  of  a  and  p  behind  the  wave 
by  the  asterisk,  *,  eq  IV.B.18.  IV.B.21  can  be  put  in  the  following  respective 

forms : 


i 


I 


The  inertial  term  given  in  eg  3  is  essentially  the  acceleration  potential  term 


P  (ty.  -  )  defined  in  eq  IX.B.44,  except  that  a  and  b  have  been  eliminated  and 

ro  O  a 

the  equation  has  been  algebraicly  reorganized.  Some  simplification  of  the 
acceleration  term  can  be  made.  Except  under  extremely  severe  shocks,  the  term 
containing  a  can  be  neglected.  Later  in  the  paper  this  will  be  done. 

In  a  steady  wave,  the  solution  to  the  balance  equations  may  be  expressed 
in  terms  of  the  single  variable 

tl  *  X-Vt.  (IV. C. 6) 


Substituting  the  constitutive  law  into  eq  IV.B.20  and  changing  variables 
from  X  and  t  to  n  by  use  of  eq  6  we  find  after  some  algebra 
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(2<VC>*C  1-  (  $Sr)) 


v2  T2  d  "  1 
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The  prime  denotes  differentiation  with  respect  to  the  variable  r| . 
integrating  the  above,  we  get 

-v2 


—  (a-o0)  (a2  -  -j  (a-K xjj 

(fi(a,a,a))  +  r  J  (  ^ 

V  I  Ja0  3  “  o-l 
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(IV.  c.  7) 
After 


(IV. C. 8) 


The  above  equation  describes  the  density  ratio  jump  across  the  wave,  once  the  wave 
speed  is  known.  However,  this  is  an  integro-differential  equation  and  would  have 
to  be  solved  by  a  numerical  method  or  some  other  approximate  method. 

In  cases  where  strain  rates  are  not  extremely  large,  eq  8  would  have 
to  be  used,  since  it  makes  use  of  the  more  complete  constitutive  equation, 
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eq  3.  Such  would  be  the  case  for  low-frequency,  large-amplitude  waves. 

However,  if  one  is  considering  shockwaves,  the  strain  rates  become 
very  large,  say  on  the  order  of  10*  -10^s  \  Under  these  extremely  high 
rates,  the  constitutive  equation  can  be  simplified  somewhat,  since  the  yield 
stress  of  ice  becomes  only  weakly  rate  dependent.  This  can  be  seen  in  eq. 

II. A. 3  where  at  large  values  of  D,  Y  changes  only  slightly  with  significant 
changes  in  D.  Consequently,  if  Y  is  replaced  by  a  constant  value  Y  =  300  bar, 

i  ° 

; 

the  constitutive  equation  developed  in  section  II. A  can  be  reduced  to  the 
following  form 


,  ,  a  ,-(Ja/a  t  d  Q(a,a,fi) 

ln  (  )e  -  2~  da 


(IV. C. 9) 


This  equation  is  obviously  more  tractable  than  the  earlier  equation.  The  value 

of  Y  =  300  bar  was  chosen  since  this  yield  stress  is  consistent  with  strains 
o 

>\,  5-1 

actually  generated  in  many  shockwaves  (D  ^  2x10  s  ) ,  and  a  range  of  D  between 

10*s_1  and  106s-1  places  the  actual  yield  stress  Y  within  ten  percent  of  the 

vA&ue  of  Y  =  300  bar. 
o 

Substituting  eq.  9  into  eq  IV.B.20,  we  obtain 
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or,  upon  integrating, 
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(IV.C.ll) 


In  steady  waves,  0=0,  and  0=0  both  ahead  of  and  behind  the  wave; 

N  • 

therefore  [Q(a,a,a)]  vanishes.  Substituting  eg.  IV.B.18  into  eg.  11  for 


steady  waves  then  gives 

a 


»*  -(2J*o /  ln  ( srr) 


-6  a/a. 
a  \  y  0 


(a*2  -  -§-  (o*-kx0))_1 

Carrying  out  the  indicated  integration  results  in 


(IV. C. 12) 


where 


p*  =  (a*  - 
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Eguations  13  and  14  give  the  pressure  jump  across  the  wave.  Eguation  10 
can  then  be  used  to  calculate  the  wave  profile,  i.e.,  the  variation  of  the  density 
ratio  a  inside  the  stress  wave.  This  can  be  accomplished  by  inverting  eg.  10  to 


obtain 


„  f  « 
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(IV. C.  15) 


where  n  is  the  variable  defined  in  eg  6,  and 
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Napadensky  (1964)  conducted  an  extensive  testing  program  on  shock  waves  in 

snow  at  Camp  Century.  These  experiments  involved  snow  with  densities  generally 

-3 

higher  than  400  kg  m  .  The  stress  waves  were  generated  by  using  a  low-density 
explosive  to  drive  a  metal  plate  into  solid  cylindrical  specimens.  By 
varying  the  amount  of  explosive  and  the  mass  of  the  driver  plate,  a  wide 
range  of  load  conditions  was  achieved.  The  explosive  technique,  however, 
usually  has  a  considerable  amount  of  experimental  scatter. 

A  streak  camera  was  used  to  record  the  motion  of  etch  marks  on  the 
snow  specimen  and  the  driver  plate.  In  this  way,  both  particle  velocity 
V  and  wave  velocity  V  were  recorded  photographically.  Since  the  speed 
of  the  camera  was  such  that  the  exposure  time  was  only  about  5  ys, 
fairly  good  time  resolution  was  achieved.  However,  as  indicated  by 
Mellor  (1977) ,  large  errors  were  incurred  in  the  experimental  program. 

There  was  a  very  poor  resolution  at  the  lower  pressure  levels,  and  it  is 
questionable  whether  or  not  Napadensky' s  reduced  data  accurately  reflected 
the  actual  stress  wave  parameters.  However,  for  the  intermediate  range 
of  stress  wave  pressures,  the  data  appeared  reasonable,  since  the  streak 
film  recordings  appeared  to  have  both  good  time  resolution  and  ability 
to  record  particle  displacement. 

Figures  33  and  34  give  a  direct  comparison  of  the  theory  and  Napadensky ’s 
experimental  results,  and  as  can  be  Been  the  two  compare  fairly  well. 
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Figure  33.  Comparison  of  theory  with  stress  wave  data 
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Figure  33  shows  the  pressure  jump  across  the  wave  as  a  function  of  the 
density  jump  Ap  caused  by  the  wave.  Figure  34  shows  the  variation  of 
plastic  wave  velocity  V  with  particle  velocity  V*.  For  particle  velocities 
above  20  n  s  \  the  data  and  theory  compare  nicely,  but  below  that,  they 
deviate  significantly.  The  theoretical  curve  terminates  at  the  point 
where  the  wave  ceases  to  be  fully  plastic. 

The  data  in  figure  34  at  the  lower  particle  velocities  are  questionable. 
Napadensky  shows  the  plastic  wave  velocity  decreasing  to  zero  as  the 
wave  intensity  (as  indicated  by  V*)  the  particle  velocity  decreases  and 
this  result  defies  physical  reasoning.  For  instance,  as  the  intensity 
of  a  shock  wave  decreases,  the  severity  of  the  plastic  deformation 
should  also  decrease.  As  the  amount  of  the  plastic  deformation  becomes 
less  significant,  the  pressure  wave  should  then  begin  to  acquire  character¬ 
istics  of  an  elastic  wave.  Equation  IV.B.20  can  be  used  to  find  the 
wave  speed: 


„  a 


-Je* 

Aa* 


(XV.C.17) 


If  Ap*/Aa*  does  not  approach  zero  as  the  stress  wave  intensity  decreases, 

2 

V  must  remain  finite  at  low  plastic  strains.  One  would  expect  Ap*/Aa*  to  acquire 
a  value  close  to  the  elastic  modulus.  This  argument  is  also  supported 
by  other  analytical  work,  such  as  that  of  Coleman  et  al.  (1964) ,  on  the 
theory  of  wave  propagation  in  nonlinear  materials.  In  particular,  we 
note  the  familiar  equation 


V2  -  E  /p.  (IV.C.18) 

P  0 

where  for  shock  waves  E  is  the  instantaneous  secant  modulus,  and  for 

r 

acceleration  s  E  -  the  instantaneous  tangent  modulus.  In  either 
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stress  wave,  E^  increases  as  the  degree  of  plastic  deformation  decreases, 
i.e.,  as  a  smaller  percentage  of  the  material  is  deformed  plastically.  As 
a  result,  appraoches  E,  Young's  modulus  of  the  material,  and  V  thereby 
increases  to  a  value  close  to  the  elastic  wave  speed  for  weak  plastic 
waves. 

In  the  case  of  strong  stress  waves  that  produce  complete  plastic 

deformation  of  the  matrix  material  and  large  jumps  in  density,  the  work- 

j 

hardening  characteristics  of  the  material  become  important  in  determining 

the  resulting  wave  speed.  In  particular,  E  increases  as  a  result  of 

P 

densification  and  work-hardening  effects,  thereby  resulting  in  larger 
balues  of  V  as  indicated  by  eq  18.  Consequently,  a  very  strong  plastic 
wave  has  a  large  wave  speed.  The  above  reasoning  justifies  the  minimum 
calculated  value  of  V  shown  by  the  curve  in  Figure  34.  This  minimum  is 
located  near  the  transition  zone  between  partially  plastic  and  fully  plastic 
waves. 

Figure  35  -  37  present  results  of  parametric  studies  of  stress  waves 
in  snow.  Figure  35  shows  how  the  density  ratio  varies  through  the  front  of 
the  stress  wave.  For  the  less  intensive  wave,  a  shorter  wavelength  is 
required  to  generate  the  acceleration  forces  necessary  to  produce  the 
compaction  to  decrease  ct  from  3.0  to  2.0. 

Figure  36  illustrates  the  variation  of  pressure  jump  (p]  with  density 
jump  for  a  range  of  initial  densities.  Note  that  as  the  initial  density  is 
increased,  the  curves  are  translated  upward.  Figure  37  also  gives  a  very 
dramatic  demonstration  of  the  effects  of  work-hardening  characteristics  of 
snow.  At  low  pressures,  the  plastic  wave  speeds  are  as  expected?  i.e., 
wave  speed  increases  with  density.  But  at  higher  pressures,  this  relationship 
is  inverted,  since  large  changes  in  a  must  accompany  the  large  pressure 


ation  of  pressure  jump  with  density  jump  and  initial  density. 


jumps  for  snow  with  initial  low  densities.  Consequently,  there  is  con¬ 
siderable  work -hardening,  which  has  a  stiffening  effect,  and  this  results 
in  increased  wave  speeds.  Snow  with  high  initial  densities  would  not 
undergo  as  much  work-hardening. 

-  Case  2;  Steady  waves  in  low  density  snow. 

For  this  case,  the  neck  growth  model  constitutive  equation  is  used, 
since  densities  as  low  as  100  Kg  m  3  may  be  studied.  Other  than  the 
constitutive  law,  the  analysis  is  identical  to  that  presented  in  Case  1. 
Consequently  only  the  results  and  discussion  is  presented  here.  The  final 
expression  for  the  pressure  jump  is 

a  _  _  (a  -^o)J 

Figures  38-40  illustrate  the  results  of  this  study.  Unfortunately 
only  a  small  amount  of  experimental  data  is  available  for  a  comparison  with 
these  results.  To  date  the  most  comprehensive  experimental  study  is  that  of 
Napadensky  (1964).  However,  Napadensky’s  experiments  were  conducted  with 
snow  with  desnities  of  0.5  Mg  m  3  or  higher,  so  a  direct  comparison  with 
the  calculated  results  here  is  not  possible  for  low-density  snow. 

Figure  38  illustrates  the  variation  of  pressure  jump  with  density  for 
several  different  initial  densities,  including  0.5  Mg  m  3,  which  is  compared 
with  Napadensky ’s  data.  As  can  be  seen,  the  results  agree  with  the  data 
for  that  initial  density.  In  Figures  38  and  39,  the  dominant  frequency  of  the 
wave  was  200  kHz.  This  frequency  would  correspond  to  what  is  produced  by 
high-speed  explosives. 

Figure  39  shows  the  manner  in  which  the  wave  speed  varies  with  the 
wave  intensity  as  determined  by  the  pressure  jump.  In  particular,  for 
each  initial  density  the  wave  speed  has  a  local  minimum  at  low  values  of 
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Figure  38.  Dependence  of  pressure  jump  on  density  jump  for  stress 
waves  in  snow  represented  by  neck  growth  model. 


wave  pressure.  This  is  attributable  to  the  work-hardening  properties  of 
snow  under  volumetric  deformations.  As  the  snow  is  compacted!  the  necks 
rapidly  thicken,  whi^h  in  part  enhances  the  material  strength.  In  addition, 
the  large  deformations  which  occurs  in  these  necks  also  work-harden  the 
neck  material,  thereby  further  increasing  the  snow  strength.  This  resulting 
increase  in  material  strength  results  in  an  increased  wave  speed  since  the 

wave  speed  is  determined  by  the  ratio  of  p*/(o*  -  <0. 

I  u 

i 

Figure  40  illustrates  the  effect  of  wave  frequency  on  the  pressure 
jump  for  pressure  waves  which  produce  a  density  jump  of  0.2  Mg  m  3.  Curves 

-3 

are  presented  for  initial  densities  of  0.1,  0.2,  and  0.3  Mg  m  .  As  can  be 
seen,  the  pressure  jump  increases  with  wave  frequency,  although  the  dependence 
of  pressure  jump  on  frequency  is  not  as  significant  as  the  dependence  on 
density. 

The  stress-wave  analysis  reported  here  has  been  shown  to  agree  well 
with  the  only  experimental  data  the  author  is  aware  of.  Unfortunately  the 
data  are  available  only  for  snow  with  an  initial  density  of  about  0.5.  Mg 
m-3.  Consequently,  no  definite  conclusions  can  be  made  about  the  validity 
of  the  theoretical  results  for  snow  with  initial  densities  below  0.5  Mg  m 
3.  However,  the  volumetric  equation  upon  which  this  study  is  based  has 
been  shown  to  represent  accurately  the  response  of  snow  under  quasi-static 
load  conditions  for  densities  as  low  as  0.1  Mg  m  3.  Consequently,  the  low- 
density  stress-wave  results  shown  here  should  at  least  be  qualitatively 
correct . 

The  wave  speed  was  found  to  be  strongly  dependent  on  initial  density 
and  the  pressure  jump  across  th  wave  front.  Since  the  pressure  is  only 
moderately  affected  by  frequency,  the  wave  speed  is  also  only  moderately 
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P  -po~0.2Mg  m~3 


pe  -0.3  Mg  m 


frequency  dependent. 

The  dependence  of  wave  speed  and  pressure  jump  on  density  is  not  as 
straight-forward  as  it  is  for  pressure  waves  which  produce  infinitesimal 
strains.  In  the  case  of  pressure  waves  that  produce  large  density  changes, 

pressure  jump  and  wave  speed  depend  on  both  the  initial  density  and  the 

* 

final  density  p  .  The  cocnplicatated  relationship  can  readily  be  observed 
in  Figure  39. 

The  type  of  pressure  wave  described  in  this  case  is  not  commonly  found 
in  nature.  However,  it  does  play  a  central  role  in  determining  the  effect¬ 
iveness  of  explosives  in  initiating  avalanches,  since  a  good  portion  of  the 
explosive  energy  that  is  transmitted  to  the  snowpack  is  absorbed  in  the 
crater  zone  through  inelastic  compaction  of  the  snow.  This  type  of  deformation 
also  is  largely  responsible  for  the  attenuation  of  the  pressure  wave  as  it 
propagates  away  from  the  crater  formed  by  the  explosive.  All  of  this 
decreases  the  energy  that  is  delivered  to  the  snowpack  for  the  purpose  of 
starting  avalanches. 

The  result  that  the  wave  pressure- jump  increases  with  wave  frequency 
would  imply  that  the  energy  absorbed  irreversibly  by  the  snow  also  in¬ 
creases  with  frequency  assuming  that  Pq  and  p*  are  held  constant.  This  in 
turn  would  imply  that  the  attenuation  rate  would  increase  with  frequency, 
since  larger  amounts  of  energy  are  dissipated  in  the  snow  as  frequency 
increases.  However,  since  this  is  weakly  frequency  dependent,  the  effect 
of  frequency  on  attenuation  rate  is  probably  not  as  important  as  p*  and  Pq 
for  plastic  waves. 
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IV. D  NONSTEADY  SHOCKWAVES 


If  the  wave  is  nonsteady,  the  wave  front  profile,  amplitude  and  wave 
speed  may  all  change  with  time.  Material  nonlinearity  and  internal  dis¬ 
sipation  mechanisms  are  accountable  for  this.  In  some  cases,  the  wave 
amplitude  may  actually  grow,  although  such  a  situation  is  generally  short¬ 
lived.  Coleman  et  al.  (1961',  have  studied  the  properties  of  nonlinear 
waves  in  some  detail.  Since  that  time  considerable  effort  has  been  devoted 
to  the  study  of  stress  wir»  in  nonlinear  or  inelastic  materials.  More 
recently,  Nunziato  afki  Walsh  (1978)  have  investigated  the  propagation  of 
waves  in  uniformly  distributed  granular  materials.  In  their  paper  they 
indicated  that,  in  a  granular  material,  the  only  density  change  induced  by 
a  shock  wave  must  come  from  compaction  of  the  matrix  material  and  not  from 
reduction  of  void  volume.  This  result  contradicts  the  results  of  this  study 
(IV.B.19)  and  the  experimental  results  of  Napadensky  (1964). 


We  have  developed  now  the  wave  equation  for  a  material  with  a  volumetric 
constitutive  equation  of  the  form  given  by  eq  IV. C. 3. 

Differentiating  eq  IV. B. 9  with  respect  to  time,  and  then  substituting 


the  equation  of  motion  eq  IV. B. 7,  yield 

i  a  /  -aP  \  _i  afa 
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(IV. D. 1) 


Then,  substituting  eq  IV. C. 3  directly  in  the  above  equation,  we  get,  after 
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(IV. D. 2) 


In  order  to  solve  this  equation,  a  finite  difference  solution  can  be  used. 


This  equation  is  extremely  nonlinear,  since  E^,  E^,  and  E2  are  all  functions 
•  •• 

of  a,  a,  and  a.  Therefore,  convergence  and  stability  problems  are  difficult 


to  handle.  A  central  differencing  technique  is  used.  At  position  and 

time  t .  let  a  have  the  value  a? .  Then  the  difference  forms  of  the  time 
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and  spacial  derivatives  of  a  become 
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The  form  shovm  in  eq  6  is  used  for  the  second  order  derivative,  since  this 
averaged  derivative  results  in  improved  stability  characteristics  of  the 
finite  difference  solution,  as  indicated  by  Ames  (1965) . 

In  solving  the  problem,  a  double  modulus  was  employed.  During  pressure 
buildup,  IV. B. 1  gives  the  appropriate  tangent  modulus,  which  becomes 
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However,  once  a  maximum  pressure  is  reached,  and  unloading  begins,  the 
static  part  of  must  decrease  in  order  to  avoid  a  large  elastic  rebound 

in  the  finite  difference  solution.  This  rebound  produces  an  artificial 

/ 

oscillation  in  the  solution  and  is  a  common  problem  with  finite  difference 
methods.  Various  methods  have  been  used  to  reduce  this  form  of  instability. 

The  easiest  way  to  facilitate  this  is  to  decrease  f^  by  a  factor  of  1/2 
when  a  becomes  negative.  A  factor  of  1/2  is  somewhat  arbitrary,  but  this 
value  was  found  to  be  sufficient  to  avoid  any  significant  volumetric  rebound. 

In  order  to  demonstrate  the  solution,  we  consider  the  particular 
problem  of  an  air  blast  directly  over  a  snow  covered  ground  surface.  We 
assume  the  blast  produces  an  overburden  pressure  of  the  form 

p  =  y  p*  (1-cos  Wt)^  0  <  Wt  <  2ir 

p  =  0  Wt  >  2TT  (IV.D.13) 

and  calculate  the  attenuation  of  the  wave  as  it  propagates  into  the  snowpack. 

The  frequency  u)  =  W/2tt  is  determined  by  the  speed  of  the  explosion,  and  the 
overburden  pressure  p*  is  determined  by  explosive  size  and  proximity  to  the 
snowpack  surface.  Of  particular  interest  is  the  response  of  a  snowpack  to 
an  overburden  pressure  of  about  20  bar  since  this  is  what  the  SULFAE  (surface- 
launched  unit  fuel  air  explosive)  weapon  system  generates. 

Figures  41-46  show  results  of  the  finite  difference  solution  of  the 
nonsteady  wave  problem  just  described.  Figures  41  and  42  illustrate  how 
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Figure  41.  Temporal  variation  of  snow  density  due  to  shockwave  loading. 


i 


i 

3 


116 


DENSITY  CHANGE  (kgra'3) 


PRESSURE  JUMP  (bar) 


Figure  44.  Attenuation  of  shockwaves  in  medium  density  snow 
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the  wave  form  is  altered  as  it  propagates  into  snow  with  an  initial  density 
of  350  kg  m-3.  The  surface  loading  has  a  frequency  of  u>  ■  5000  Hz  and  an 
amplitude  of  p*  =  12  bar.  Figure  41  shows  the  temporal  variation  of  the 
density  change  for  points  at  0,  1,  2,  5,  and  8  cm  into  the  snowpack. 

Figure  42  demonstrates  the  density  profile  for  various  times.  Note  how  the 
wave  spreads  and  the  density  rates  p  decreases  as  the  wave  propagates 

further  into  the  snowpack. 

i 

Figures  43  and  44  make  a  direct  comparison  of  three  different  pressures. 
Figure  43  compares  wave  attentuation  in  terms  of  density  jump  and  Figure  44 
shows  the  pressure  attentuation.  As  can  readily  be  seen,  the  advantage  of 
the  higher  pressures  is  largely  eliminated  within  the  first  10  cm.  This 
merely  points  out  the  substantial  energy  absorbing  capability  of  snow. 

As  expected,  the  highly  dissipe'  Lve  characteristics  of  snow  rapidly 
change  the  stress  wave  as  it  propagates  through  snow.  For  pressure  waves 
with  a  magnitude  in  excess  of  5  bar,  the  pressure  amplitude  reduced  to  a 
small  fraction  within  10  cm.  This  result  is  to  some  extend  verified  by 
Wisotski  and  Snyder  (1966) .  In  the  tests  reported  by  Wisotski  and  Snyder, 
one-pound  spherical  Penolite  charges  were  detonated  in  deep  mid-season 
snowpack.  Piezoelectric  gages  were  used  to  record  arrival  times  and  pressures 
to  within  0.15  m  of  the  charge.  These  transducers  apparently  had  a  broad, 
flat  frequency  range,  so  that  the  recorded  results  should  be  meaningful. 

They  also  observed  that  the  wave  speed  close  to  the  charge  was  significantly 
slower  than  further  from  the  charge.  However,  there  was  a  great  deal  of 
scatter  in  measurement  of  arrival  times,  so  no  precise  measurements  were 
made.  But  the  pressure  readings  for  a  6-bar  pressure  wave  showed  good 
agreement  with  the  results  shown  in  Figure  44. 

Figure  45  shows  the  variation  of  the  stress  wave  speed  as  the  wave 
propagates  into  the  snow.  Initially  the  21-bar  wave  travels  at  a  signi- 
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ficantly  higher  speed,  but  this  situation  quickly  changes.  The  variation 
in  stress  wave  speed  is  due  to  a  combination  of  factors.  The  wave  speed 
initially  decrease  as  the  waves  attenuate  and  work-hardening  effects  become 
less  significant.  However,  once  the  wave  intensity  is  reduced  to  a  critical 
value,  dissipative  effects  likewise  become  less  significant  and  the  wave 
speed  begins  to  increase.  These  results  are  in  close  agreement  with  those 
shown  in  Figure  37  for  steady  waves. 

I 

Figure  46  shows  the  effect  of  wave  frequency  on  attenuation.  One  can 
readily  see  that  higher  frequency  waves  do  attenuate  more  quickly,  but  this 
effect  is  not  all  that  great.  Of  the  three  waves  shown,  the  difference  in 
pressure  jump  by  the  5_cm  position  is  only  about  20%.  The  frequency  difference 
is  partly  erased  by  wave  spreading,  and  probably  the  higher  frequency  wave 
spread  more  quickly  and  therefore  experience  a  quicker  reduction  in  frequency 
content'. 

It  is  also  possible  to  study  the  growth  and  decay  of  shock  waves  by 
considering  the  jump  equations.  These  equations  cannot  give  as  much  information 
as  the  direct  approach  just  discussed,  since  stress  wave  profile  and  wave 
length  cannot  be  calculated  from  the  jump  equations.  However,  the  direct 
approach  using  the  finite-difference  technique  is  computationally  a  time- 
consuming  and  expensive  means  of  solving  the  problem.  If  wave  attenuation 
is  sought,  there  should  be  more  convenient  means  of  doing  this. 

IV. F.  THE  USE  OF  JUMP  EQUATIONS  FOR  ANALYZING  NONSTEADY  WAVE  PROPAGATION. 

A  solution  to  the  problem  of  shockwave  propagation  in  a  nonlinear 
dissipative  material  can  be  a  complicated  problem  if  one  attempts  to  determine 
such  wave  properties  as  wavespeed,  pressure  profile,  attenuation,  and 
alteration  of  the  wave  as  the  wave  propagates  into  a  medium.  Spence  (1973) 
characterized  wave  propagation  properties  for  viscoelastic  materials. 
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Generally  a  detailed  analysis  requires  solution  of  the  momentum  and  mass 


balance  equations.  In  particular,  if  one  desires  a  quantitative  evaluation, 
a  numerical  solution  of  these  equations  is  usually  required,  but  this  can 


be  an  expensive  and  time  consuming  procedure,  since  the  solution  of  the 
differential  wave  equation  can  involve  a  large  number  of  simultaneous 
difference  equations. 

This  can  best  be  illustrated  by  considering  the  wave  equation  solved 


numerically  in  the  preceding  section: 


_  32q  .  3ET  3a  o  32a  .  3  3a  ^  „  3ax,T„„  ,4 

ET  a  2  3X  3X  “  a  ,  2  3X  (E1  3x  E2  3x^ {IV£  *1) 

3x  0  3t 

The  constitutive  law  (for  instance  eq  IV. C. 3)  can  be  conveniently  written 


in  the  following  form 


p  -  Q, (a,aQ)  +  Q2(a,aQ)a  +  Q3(a,a  Q)a 


(IV.  Ji2) 


The  functions  Q^,  Q2  and  can  readily  be  determined  from  eq  IV. C. 3.  For 

purposes  of  this  discussion,  it  is  convenient  to  use  eq  2.  E^,  E^,  and  E2 

•  .. 

are  functions  of  a, a, a,  so  that  the  wave  equation,  which  is  basically  the 
differential  equations  of  motion,  and  continuity,  is  a  nonlinear  fourth 
order  partial  differential  equation,  g  ,  and  are  all  continuous  and 
differentiable  functions  of  a  and  the  initial  ratio,  aQ.  The  above  equation 
was  shown  to  accurately  characterize  the  pressure  response  of  snow  to  high 
rate  compressive  deformations. 

In  this  section  a  solution  to  the  stress  wave  problem  is  presented  and 
compared  to  earlier  results.  The  method  used  here  entails  the  use  of  jump 
equations,  often  referred  to  as  weak  solutions.  The  shockwave  is  defined 
as  a  smooth  one  parameter  family  of  points  Y(t)  such  that  Y(t)  gives  the 
material  point  at  which  the  stress  wave  is  located  at  time  t.  X  is  the 
position  of  the  particle  in  the  reference  configuration,  and  x  ■  x(X,t)  is 
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the  particle  position  for  any  time  t.  The  intrinsic  velocity  of  the  wave 
is 


V  "  It  Y(t)  "  *  (IV.  8.3)  £ 

If  f  (X,t)  is  any  descriptive  variable  of  the  material,  the  jump  in  f, 
denoted  as  [f]  is  the  difference  in  the  values  of  f  just  behind  and  in 
front  of  the  wave: 

If)  -  t  -  f+  (IV.  £.4)  £ 

where 

f-  *  lim  f (X-6,  t) 

6-*  o 

f+  «  lim  f  (X+6,t) 

6  -*■  o 

Here  we  assume  the  wave  is  propatating  in  the  positive  X  direction.  This 
characterization  of  a  shockwave  depicts  the  wave  as  a  singular  surface, 
i.e.  the  wave  is  assumed  to  have  an  insignificant  wavelength,  so  that  the 
only  properties  which  can  be  garnered  from  this  description  are  the  jumps 
in  variables  across  the  wave.  As  a  consequence,  some  information  is  lost, 
since  the  internal  structure  of  the  stress  wave  cannot  be  studied.  However, 
the  problem  becomes  mathematically  much  more  manageable. 

In  the  following  the  governing  equations  and  widely  recognized  relations 
for  stresswaves  are  again  presented  for  more  convenient  reference.  Then 
additional  results  are  developed  to  make  possible  calculations  of  the 
growth  and  decay  of  stress  waves  by  use  of  the  jump  equations.  Finally, 
the  current  development  is  compared  to  earlier  results. 
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In  what  follows ,  we  assume  the  stress  wave  is  propagating  into  an 


undisturbed  medium.  Thus  a  *  <xq,  and  the  pressure  p,  particle  velocity 
v,  and  density  ratio  rate  a  ahead  of  the  wave  all  vanish.  Therefore, 


la]  -  a  -a 


(XV. B. 5) 


Ipl  “  P 
tv]  =  v 


[a]  -  a 


The  derivative  of  the  pressure  p  with  respect  to  either  time  or  position 
has  the  form 


-E  a  +  E,  a  +  E„  a 
T  1  2 


P  -  -e  a  +  E,a  +  E_a  , 
x  T  x  lx  2  x 


(IV.  E6) 


where  a  subscripted  x  denotes  differentiation  with  respect  to  X.  Then 
the  jumps  in  these  variables  may  be  shown  to  be 


[p]  -  -  e"  [a]  +  e~  [a]  +  e"  [*a*J 


(IV.E  .7) 


[Pxl  -  -  et  taxl  ♦  laxl  +  b2  tax] 

ET,  E^,  and  E2  are  the  values  of  E^,  and  just  behind  the  wave.  The 
rate  of  change  of  [p]  with  respect  to  time  may  also  be  shown  to  be 


cpi t  =  i«it  +  Ei  i«jt  +  e2  i5it 


(IV. ^.8) 


The  subscripted  t  implies  differentiation  with  respect  to  time.  In  the  above 
Ip] represents  the  change  or  jump  in  the  value  of  the  pressure  derivative  p 
across  the  wave,  whereas  £p J gives  the  rate  of  change  of  the  pressure  jump 
(pi  as  the  wave  propagates  through  the  medium. 

For  a  plane  wave  propagating  in  the  X  direction,  the  two  governing  differ¬ 
ential  equations  are  the  momentum  and  mass  balance  equations.  These  are; 


JE.-c 


(IV.E. 9) 
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3v 

3x 


/a 


(IV.E.10) 


these  equations  may  be  put  in  the  form  of  jump  equations  in  much  the  same 
manner  as  was  done  by  Coleman  (1964)  and  Chen  and  Gurtin  (1971) .  The  difference 
form  for  these  equations  are: 


I  3p/3x)  “-(»otv] 


I3v/3x] *  Ia/ao] 


Alternate  forms  of  these  two  equations  are 


I  Pi  -  PQV[v] 


lv]  *=  -V[  3x/3x]  =  -  V]F] 


I  v]  -  -  [a] 

o 


F  is  the  deformation  gradient,  3x/3x. 

Equations  13  and  15  may  be  used  to  arrive  at 

v2  .  _!°JeI 
pq  la] 

The  secant  modulus  for  the  material  is  defined  as 


e"  .  -  SeL 

a  (a] 

so  that  the  wave  speed  V  then  has  the  form 


V  -  o  E/p 
o  s  o 


(IV.E.il) 

(1V.E.12) 


(IV.E.13) 

(1V.E.14) 

(IV.E.15) 


(IV.E-16) 


(IV.E.17) 


(IV.E.18) 

Zn  addition  to  the  above,  the  familiar  compatibility  relation  for  shockwaves 
(Coleman  et  al  1964)  is  of  considerable  use  in  the  study  of  shockwaves.  If 
f(X,t)  is  any  function  of  position  and  time,  the  compatability  relation  requires 
that 

[f)t  -  If)  +  VIfxl  (IV.E.19) 
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The  above  relations,  eqs  11  ~  19,  are  considered  to  be  the  standard  jump 
equations  for  the  propagation  of  shockwaves  in  materials.  They  have  been  used 
extensively  in  previous  studies  by  Coleman  et  al  (1964) ,  Carrol  et  al  (1973) , 

Chen  et  al  (1971) ,  and  others  to  characterize  the  manner  in  which  shockwaves 
are  propagated  through  materials,  irrespective  of  the  material  properties. 

We  now  use  the  foregoing  relations  to  evaluate  the  growth  and  decay  rates 
of  shockwaves  in  snow.  Differentiating  eq  15  and  using  eq  19  with  f  =  v 
results  with  l 

Via]  +710],.  *  -a  dv]  +  v[v  ])  (iv,e.20) 

to  X 

Then  substituting  eqs  11  and  12  gives 

a 

V  la]  +  V[a]  =  r2-  [p  ]  -  V[a]  (IV.E.21) 

t  pQ  x 

This  equation  is  then  further  modified  by  using  eq  7,  eq  9  with  f  =  a  and 
f  «*  a,  and  eqs  13  and  15.  After  some  algebra,  we  get 

o  PV  tcxl  .  _  . 

- -  V  »  -  (ET  +  Es)  la]t  +  (Et  -  Eg)  [a]  (IV.E.22) 

o 

+  E~  ([a]t  -  [a])  +  E~  (-fa)t  -  fa* J > 

V  may  be  eliminated  by  differentiating  eq  15  and  substituting  into  eq  22  for 
V.  The  result  is 

|(3e;  +  E~)  [alt  -  ~  E1ta]t  =  (e-t  -  e~) [a] 

-  E1la]  +  e"  (  |  fa] t  -  [a*])  (iv.e.23) 

As  can  be  seen  here,  the  full  collection  of  moduli  determine  the  rate  of 

l 

change  of  la]  for  the  above  form  of  the  wave  equation. 

Some  additional  criterion  for  determining  what  constitutes  the  wave 
amplitude  is  needed.  For  instance,  one  can  require  that  the  wave  has  passed  a 
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point  where  the  material  is  no  longer  being  compressed,  i.e.  when  a®  o,  which 
would  be  equivalent  to  requiring  [a]  =  o.  Rather  than  this,  we  impose  the  re¬ 
quirement  that  the  wave  has  passed  a  point  when  the  peak  pressure  has  been 
reached.  This  is  equivalent  to  the  condition 

[p]  -  o  (IV  J3  .24) 

This  then  implies,  with  eq  7 

-e"  [a]  +  e"  [a]  +  E~  (*a*l  *■  o  (iv.E.25) 

With  this,  [a]  can  be  eliminated  from  eq  23,  and  we  get 


i(3E~  +  e~) [a?t  -  ~  lo*t  *  -e"  [a)  +  j  e"  [a)t  (iv.e.26) 

Eq  26,  in  conjunction  with  eqs  2,  and  8  may  be  used  to  study  the  manner  in 
which  the  stress  wave  propagates  through  snow.  Given  an  initial  set  of  conditions, 
say  [a]  and  (a],  these  equations  may  be  solved  to  determine  [p]fc,  [aJtand 
IO)t  and  by  numerical  integration,  the  temporal  variation  of  [p] ,  [a] , 
and  (a]  can  be  evaluated. 

Results  in  Section  II  indicate  have  the  acceleration  modulus  is  extremely 
small,  and  the  terms  E 2  [a)  in  eqe  8  and  26  may  be  neglected.  Therefore  we 
will  assume  that  effects  associated  with  E 2  can  be  neglected  and  drop  this 
term  from  eq  26. 

The  formulation  developed  here  may  now  be  compared  to  results  obtained  in 
the  earlier  development  where  Eq  1  was  integrated  by  a  finite  difference 
method.  In  that  section,  the  material  constitutive  equation  for  volumetric 
deformation  was  given  by  eq  2  where  Q^,  Q2,  and  Q3  are 

2Yo  j 

Qx  *  — a  ■  exp(-<pa/aQ)  In  (a/(a-l) 


C2 


(a-1) 


-4/3  _a-4/3} 


(IV.E.27) 


03  -  -  ~  (  (a-l)"1/3  mx‘1/3) 


Xn  the  earlier  analysis,  it  was  determined  that  a  finite  difference  solution 


to  the  wave  equation  gave  very  reasonable  results  when  compared  to  what  experi- 

\vJ 

mental  results  are  available.  However,  as  is  common  in  finite  difference 
solutions  of  higher  order  differential  equations,  stability  problems  were 
incurred  and  had  to  be  solved,  and  the  required  grid  size  resulted  in  long 
computation  terms. 

Eq  26,  with  the  term  with  eliminated,  can  readily  be  put  in  a  dif¬ 
ference  from  in  order  to  obtain  expressions  for  [a]t  and  [a] t  at  any  instant, 
t.  This  equation  and  the  compatability  equation  eq  19  with  f  =  a  can  then  be 


used  jointly  to  solve  for  [a]^  and  [a] ^  The  compatability  condition  is 


[a]t  *>  [a]  +  v  [axl 


(IV.®. 29) 


can  be  evaluated  by  expanding  a  into  a  Taylor  series  about  a  point  just 
behind  the  shockwave.  To  demonstrate  this,  consider  Figure  47,  which  shows 
two  aeprofiles  resulting  from  the  shockwave.  The  solid  line  describes  the 
wave  profile  at  a  time  ,  and  the  dashed  line  represents  the  wave  profile  at 
time  t r'  *  tJ  -At,  where  At  is  a  small  time  increment.  The  wave  has  progressed  a 
distance  Of  Ax  «  vAt  during  this  time  from  position  X^  ^  to  X..  °j  l  rePresents 
the  density  ratio  at  time  t?  and  position  X..,^-=  Xj  ~  VAt.  A  Taylor  series 
expansion  of  a  about  point  ^  is  of  the  form 


^  I —  a  (tj_1,  X.  )  (At)n 

j-1  n-o  nl  3^n  j-1 


which  may  be  truncated  to 


*3-1 -«&  ♦i3-i  At 


DENSITY  RATIO, 


X- COORDINATE 


Figure  47.  Idealized  density  profiles  in  a  shockwave. 


for  small  values  of  At.  Then,  V [a^]  can  be  shown  to  have  the  form 


(c^  -  ) 

via  ]  -  v  IV  J-r 


Ax 


aj  -  a3'1  -  aj_1  At 

-  aj  aj-l  3-1 

At 


and  eq  29  then  has  the  difference  form 


la]; 


•  1  aJ 

j  +  ^ 


aj'}  -  aV]  At 

„.l-A _ Hi _ 

At 


(IV.E.30) 


and  eq  26  in  difference  form  becomes 

i  (3(E“)^  +  (E~)  j)  R*]^  -  f(E~)^[a]^  =  -  (e;>^  !<*]  j  (IV.  E.  31) 

•  j 

Upon  using  eqs  30  and  31  to  evaluate  faJt  and  (alt  at  any  instant  t  where  the 
wave  is  located  at  position  X..,  then  [a]  ana  [a)  can  be  calculated  by  direct 
step  forward  integration. 

Figures  48  -  51  illustrates  the  results  of  the  method  described  above. 

In  figure  48,  wave  pressure  attenuation  in  medium  density  snow  is  calculated 
and  compared  to  experimental  results  of  Wisotksi  and  Snyder  (1964) .  As  can  be 
seen,  agreement  with  the  data  is  excellent,  although  one  should  notice  that 
there  is  about  a  twenty  five  percent  difference  between  the  density  used  in 
the  calculations  and  in  snow  studied  by  Wisotski  and  Snyder.  It  is  anticipated 
that  the  attentuation  in  the  lower  density  snow  should  be  larger.  Unfortunately 
there  still  is  a  shortage  of  good  shockwave  data  for  making  a  more  rigorous 
evaluation  of  the  theory. 

Figure  49  illustrates  the  attenuation  of  the  pressure  jump  due  to  the 
shockwave.  Also  included  in  the  figures  are  results  obtained  by  the  full  numerical 
solution  of  eq  1  in  the  preceding  section  and  in  the  paper  by  Brown 
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Figure  48,  Comparison  of  theoretical  solution  to  experimental  data. 
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Figure  49.  Comparison  of  two  methods  for  calculating  shockwave  attenuation 
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Figure  50.  Comparison  of  two  theoretical  solution  techniques 
(Sections  IV  C  (  D). 
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(1980b)  and  in  Section  IV. C.  As  can  be  seen,  some  discrepancies  between  the 
two  theories  emerge.  There  appears  to  be  a  fairly  good  correlation  between 
the  two  theories  for  the  lower  pressure  waves,  but  for  a  pressure  wave  with  an 
initial  amplitude  of  21  bar,  the  two  solutions  diverge  rather  quickly.  It  is 
not  yet  detennind  which  theory  is  more  correct.  While  the  direct  finite 
difference  solution  of  eq  1  would  represent  a  more  rigorous  method  than  the 
method  presented  here,  the  stability  problems  associated  with  numerical  solution 
of  high  order  difference  equations  may  introduce  some  numerical  errors. 
Unfortunately  there  is  no  experimental  data  with  which  to  compare  these  results. 

Figure  50  compares  the  attenuation  in  the  density  jump,  and  Figure  51 
illustrates  how  the  wave  speed  changes  as  the  wave  propagates  through  a  homo¬ 
geneous  snowpack.  In  both  of  these  figures  there  do  exist  differences  between 
the  two  formulations  although  the  results  do  show  very  similar  trends. 

The  method  described  here  for  characterizing  shockwave  propagation  in 
materials  such  as  snow  has  not  been  studied  extensively  in  the  past.  When 
jump  equations  have  been  used,  their  use  has  been  restricted  almost  entirely 
to  calculations  of  density  jump  or  wave  speed.  Very  little  effort  has  been 
made  toward  evaluating  wave  attentuation  or  growth. 

As  has  been  shown,  the  use  of  weak  solutions  do  appear  to  be  a  viable 
means  of  investigating  shockwaves  without  having  to  resort  to  complicated  and 
expensive  numerical  solution  of  the  wave  equation.  While  a  numerical  solution 
was  employed  here,  effort  required  to  obtain  a  solution  was  but  a  small  frac¬ 
tion  of  the  effort  required  for  the  solution  of  eq  1  in  the  paper  by  Brown 
(1980b) . 

The  numerical  results  here  do  imply  that  snow  does  strongly  attenuate  in 
plastic  stress  waves.  As  has  been  shown,  the  wave  amplitude  is  reduced  to 
about  10%  of  its  original  amplitude  within  10  cm  of  propagation. 
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V.  EXPERIMENTAL  PROGRAM  ON  STRESS  WAVES  IN  SNOW. 

Other  than  the  pioneering  work  by  Napadensky  (1964) ,  there  has  been  very 
little  significant  laboratory  work  on  plastic  shock  waves  in  snow.  Other  work 
includes  studies  on  shock  waves  in  wet  snow  by  Wakahama  and  Sato(1977) .  Brown 
(1979(a),  [b] ,  [c]  has  made  theoretical  studies  of  steady  and  non-steady  waves 
in  snow,  and  the  results  appear  reasonable  when  compared  to  what  little  experi¬ 
mental  evidence  is  available  for  stress  waves. 

In  order  to  improve  on  the  quantity  and  quality  of  shock -wave  data,  an 
electromagnetic  stress-wave  generator  has  been  constructed,  and  this  instrument 
is  described  in  detail  here.  Other  methods  of  generating  shock  waves  could 
have  been  used.  For  instance,  explosives  can  be  used  to  produce  a  wide  range 
of  wave  amplitudes,  although  wave  frequency  is  not  an  easily  controlled  parameter 
and  explosives  also  usually  have  a  great  deal  of  experimental  scatter.  The 
electromagnetic  generator  was  chosen,  since  it  appeared  to  provide  an  optimum 
combination  of  pressure  and  frequency  capability.  In  addition  it  has  negligible 
experimental  scatter  and  is  easy  to  use.  One  drawback  is  the  danger  associated 
with  high  voltages,  and  such  a  system  must  have  incorporated  into  it  safety 
features  to  mxnimize  this  aspect  of  the  system.  However,  the  use  of  explosives 
also  requires  implementation  of  safety  features. 

Basically,  the  electromagnetic  generator  consists  of  a  20  kV  capacitor 
with  a  capacitance  of  30  WF,  a  triggering  unit,  a  spark-gap  trigger,  and  a 
system  of  conducting  plates  and  load  strips.  The  generator  works  on  the 
principle  that  when  two  high-energy  electric  currents  are  conducted  in  opposite 
directions  along  adjacent  parallel  strips,  a  strong  repulsive  electromagnetic 
force  is  generated  to  force  the  two  strips  apart. 

The  capacitator  and  high-voltage  side  of  the  trigger  are  both  completely 
submerged  in  a  tank  filled  with  transformer  oil,  thereby  insulating  components 


L38 


and  reducing  any  corona  effects  that  might  occur  during  charging  of  the  capacitor. 
A  20  kV  regulated  power  supply  is  used  to  charge  the  capacitor,  and  charging 
time  generally  is  less  than  S  min.  The  capacitor  voltage  can  be  read  from 
either  a  gauge  on  the  instrument  panel  or  from  a  digital  voltmeter  if  more 
accuracy  is  required.  A  400  Mil  dropping  resistor  is  used  to  drop  the  voltage 
to  the  gauge.  The  digital  voltmeter  is  connected  to  opposite  ends  of  a  0.4  Mft 
resistor  in  series  with  the  dropping  resistor,  so  that  a  1  V  output  on  the 
digital  meter  corresponds  to  a  1  000  V  charge  in  the  capacitor.  The  digital 
voltmeter  is  connected  remotely  to  the  generator,  so  that  charge  voltage  can 
be  monitored  from  outside  the  cold  room  in  which  the  testing  is  done. 

As  can  be  seen  in  figure  52,  the  electrical  circuitry  consists  of  three 
separate  components:  (1)  the  trigger  control  system;  (2)  the  power  supply, 
capacitor,  and  conducting  elements;  and  (3)  the  monitoring  instruments.  While 
these  three  separate  systems  are  interconnected,  they  serve  different  purposes. 

The  trigger  control  unit  outlined  in  Figure  52  serves  several  purposes. 

First  it  controls  power  to  the  entire  system.  Second,  when  the  power  is 
switched  on,  an  electromagnet  activates  a  high-voltage  switch  to  close  the 
charging  circuit  to  allow  the  capacitor  to  be  charged.  If  the  system  is 
turned  off  or  a  power  failure  occurs,  the  switch  automatically  opens  to  allow 
a  limiting  resistor  to  discharge  the  capacitor.  In  addition,  a  100  Mfi  bleed 
resistor  is  also  permanently  installed  in  the  circuit  to  bleed  the  capacitor 
down  in  case  the  H.v.  switch  for  the  other  resistor  should  fail.  The  control 
unit  also  controls  the  power  to  the  high-voltage  power  supply  and  contains 
the  trigger  module  which,  by  means  of  a  transformer,  pulses  the  spark  gap  with 
a  25  kV  pulse.  This  ionizes  the  gases  in  the  trigger  which  closes  the  circuit 
and  allows  the  current  to  be  conducted  across  the  spark  gap,  through  a  control 
coil  to  the  load  strip  and  to  the  ground  side.  The  coil  is  used  to  make  the 
inductance  of  the  complete  system  large  enough  so  that  the  current  will  not 
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exceed  the  maximum  allowable  current  for  the  spark  gap. 


V. A  ELECTROMAGNETIC  THEORY 

The  generator  is  a  LRC  circuit,  and  the  system  performance  is  determined 
by  the  differential  equation 


S  “*>  *RiS*§-° 

at 
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where  Q  is  the  instantaneous  capacitor  charge,  L  is  the  circuit  inductance,  R 
is  the  total  circuit  resistance,  and  C  is  the  capactitance.  In  the  underdamped 
mode,  the  solution  to  eq  1  is 


Q(t)  *=  [Qq  exp  (  -  at) 
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sin  tut  +  cos  tut),  (V.A.2) 


a  -  R/2L, 

W  =  1/LC  -  (R/2L)2 


(V. A. 3) 
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The  instantaneous  current  i(t)  can  be  shown  to  be 

Co 

i(t)  *  exp  (-at)  sin  tut.  (V.A.5) 

The  current  is  exponentially  damped,  and  if  T  is  the  period  of  the  oscillating 
current,  the  full-cycle  peak-to-peak  ratio  of  the  current  amplitude  is 


—  =  exp  (  -  aT) 
A1 
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where  A^  and  A^  are  the  peak  current  values  of  two  successive  complete  cycles. 
The  values  of  L  and  C  may  then  be  solved  for  in  terms  of  C,  Ty  and  Aj/A^. 


This  yields  (Snell  and  others,  1973) 
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Hie  pressure  p  produced  between  the  two  conducting  strips  varies  with  the 


current  as 


p 


(V.A.9) 


where  w  is  the  strip  width,  and  u  is  the  magnetic  permeability  of  the  material 
separating  the  strips,  which,  in  this  case  is  0.254  mm  thick  Mylar.  Equation 
9  is  valid  so  long  as  the  strip  width  w  is  large  compared  to  the  separation  of 
the  conducting  plates. 

V.B  ASSESSMENT  OF  GENERATOR  PERFORMANCE 

In  its  present  form,  4he  generator  has  about  6000  J  of  energy  when  charged 
to  20  kV.  Maximum  current  with  this  configuration  is  about  160  kA.  The  spark 
gap,  however,  is  designed  for  a  maximum  current  of  about  100  kA,  so  the  maximum 
voltage  used  to  date  is  about  15  kV.  This  level  of  current  and  voltage  was 
however,  found  to  be  sufficient  to  produce  high-pressure  shock  waves  in  snow. 

When  the  system  was  originally  completed  and  first  tested,  the  resistance 
and  inductance  were  found  to  be  3.75  mfl  and  85  nH,  respectively,  and  a  capacitor 
voltage  of  15  kV  resulted  in  currents  of  about  275  kA,  which,  after  about  20 
shots,  destroyed  the  spark  gap.  After  that,  the  control  coil  shown  in  figure  52 
was  placed  in  series  between  the  capacitor  and  load  strips  to  raise  L  and  R  to 
values  of  650  nH  and  19  m^,  which  keeps  the  maximum  current  below  100  kA  when 
the  capacitor  voltage  is  less  than  15  kV.  This  is  significantly  lower  than 
the  potential  system  output  (450  kA  at  20  kV)  when  the  control  coil  is  removed, 
but  to  date  a  spark  gap  to  handle  these  current  levels  has  not  been  constructed 
for  the  stress-wave  generator. 

The  current  in  the  circuit  is  monitored  by  a  Rogowski  coil  that  is  encircled 
by  the  conducting  strips  near  the  end  of  the  load  strips.  This  coil  is  simply 
a  length  of  insulated  wire  that  is  looped  and  twisted  tightly  into  a  shorter 
length.  When  this  is  placed  adjacent  to  the  conducting  strip  (in  this  case, 
at  a  point  near  the  load  strip) ,  the  electromagnetic  field  set  up  by  the  large 
current  produces  a  voltage  across  the  leads  to  the  Rogowski  coil.  Actually  it 
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makes  little  difference  if  the  coil  encircles  the  strip  or  is  encircled  by  the 
strip  itself,  as  was  done  in  our  case.  The  electromagnetic  field  produces  a 
voltage  across  the  coil  leads  in  either  case.  The  voltage  across  the  leads 
from  the  Rogowski  coil  (Snell  and  others,  1973)  is  proportional  to  the  instan¬ 
taneous  voltage  in  the  capacitor.  The  output  from  a  typical  oscilloscope 
trace  of  the  coil  output  is  shown  in  figure  53.  Snell  and  others  (1973)  have 
shown  that  the  maximum  current  can  be  determined  from  the  period  T  of  the 
Rogowski  output  and  the  full-cycle  peak-to-peak  ratio  A^/A^  of  the  voltage 
generated  across  the  coil.  This  relation  is 


From  figure  53,  A  /A  =  0.611,  T  =  30  ys,  so  for  a  capacitance  of  30  yF 
and  a  voltage  of  15  kV,  a  maximum  current  of  84  kA  is  developed.  Under  these 
conditions,  a  pressure  in  excess  of  5  000  kPa  can  be  produced  for  a  specimen 
with  a  10  cm  x  2.5  cm  cross-section. 

The  test  procedure  is  relatively  simple.  While  the  capacitor  banks  are 
uncharged,  the  specimen  is  cut  and  shaped  and  mounted  on  the  load  strips,  and  the 

pressure  transducer  is  then  positioned  at  the  top  end  of  the  specimen.  Then  the 

capacitor  bank  is  charged,  which  generally  takes  less  than  5  min.  When  the  de¬ 
sired  voltage  is  reached,  the  spark  gap  is  pulsed,  thereby  completing  the  test*-. 
The  Rogowski  coil  and  transducer  output  are  recorded  on  trace  oscilloscopes. 
Repeatability  of  the  system  is  excellent,  with  generator  output  scatter  being 
less  than  5%. 

V.C.  TEST  RESULTS  WITH  THE  ELECTROMAGNETIC  STRESS  WAVE  GENERATOR 

A  testing  program  was  initiated  during  the  research  project  funded  by  the 
Army  Research  Office.  The  purpose  of  the  project  was  to  evaluate  the  propagation 

of  shockwaves  in  low  density  snow.  In  all,  about  150  shockwave  tests  were  run  on 
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specimens  with  initial  densities  ranging  from  150  Kg  m  to  300  kg  m  .  The 
reduced  data  is  represented  in  Figures  54-58.  These  tests  are  separated  into 
ranges  of  density  and  the  pperSting  voltage  of  the  stress  wave  generator. 

As  was  indicated  earlier,  the  spark  gap  trigger  was  incapable  of  withstanding 
currents  in  excess  of  100  kA,  even  though  the  stress  wave  generator  had  the 
capability  of  producing  currents  as  large  as  460  kA.  Therefore,  the  tests  had  to 
be  run  at  pressures  which  were  only  about  5%  of  its  full  capability.  Financial 
constraints  would  not  allow  the  purchase  of  a  trigger  which  could  handle  these 
large  currents. 

In  all  tests,  natural  snow  gathered  from  the  Bridger  Mountains  was  used. 

Seived  snow  normally  has  densities  in  excess  of  350  kg  m  ^  and  hence  could  not 
be  used.  In  addition  seived  snow  has  different  mechanical  properties  than  natural 
snow  of  comparable  densities,  since  seiving  alters  the  granular  bonding  and  the 
manner  in  which  the  grains  interact  with  each  other  during  loading  and  deformation. 
Since  interest  was  in  natural  snow,  it  was  decided  to  not  use  seived  snow. 

The  decision  to  use  natural  snowpack  resulted  with  its  own  set  of  problems. 

For  instance,  natural  snowpack  is  normally  nonhomogeneous.  Density  tends  to  vary 
considerably,  and  it  is  also  difficult  to  avoid  finding  small  ice  crusts,  surface 
hoar  layers  or  sun  crust  layers  in  the  snow.  Sometimes  these  layers  are  very  hard 
to  detect.  However,  they  do  affect  the  nature  in  which  stress  waves  are  propoa- 
gated  through  snowpack.  They  act  as  a  surface  which  reflects  part  of  the  wave 
if  the  surface  is  perpendicular  to  the  propagation  direction. 

Natural  snow  is  also  much  more  difficult  to  handle  and  test  than  seived 
snow,  particularly  if  the  density  is  low.  Often  specimen  damage  can  occur,  thereby 
sesulting  with  bad  test  results.  In  addition,  it  is  often  difficult  to  obtain 
good  surface  contact  between  the  specimen  and  loading  strip  and  betweeen  the  specimen 
and  the  pressure  transducers. 


The  difficulties  cited  above  results  with  an  abnormally  large  number 
of  bad  shockwave  tests.  Of  the  approximately  150  tests  that  were  run,  only 
a  little  more  than  50  tests  gave  valid  results,  and  these  are  shown  in 
Figures  54-58.  These  tests,  however,  do  form  a  basis  upon  which  to  evaluate 
shockwave  propagation  properties  of  low  density  snow. 

At  the  time  of  the  writing  of  the  report,  an  analytical  solution  of 
nonsteady  wave  propagation  in  low  density  snow  was  being  worked  on  but  as 
of  then,  unfinished.  This  solution  was  was  based  on  the  neck  growth  model 
described  earlier  in  this  report. 

The  test  data  does  show  that  shockwaves  do  attenuate  rapidly  in  low  density 
snow.  Shockwave  amplitude  drops  to  less  than  10  percent  within  10  cm  of  propa¬ 
gation  distance.  These  results  are  consistent  with  results  obtained  previously 
for  higher  density  snow. 


VI.  CONCLUSIONS 


Presented  in  this  report  has  been  the  research  efforts  of  the  author  during 
the  three  year  period  of  a  research  grant  funded  by  the  Army  Research  Office, 

Research  Triangle  Park,  NC.  The  purpose  of  the  grant  was  the  investigation  of 
high  rate  volumetric  properties  of  snow  with  special  emphasis  on  shockwaye  propa¬ 
gation.  It  is  felt  that  progress  was  made  in  gaining  a  better  understanding  of 
the  unique  properties  of  snow  when  subjected  to  high  rate  deformations.  Applications 
of  theory  to  such  problems  as  vehicle  mobility  in  snow  and  shookwave  propagation 
has  been  shown  to  give  results  which  seem  realistic  and  reasonable.  This  last 
phrase  was  used, since  there  still  does  not  exist  enough  experimental  data  to  draw 
definite  conclusions. 

Hore  work  still  needs  to  be  done  to  gain  a  better  understanding  of  the  high 
rate  properties  of  snow.  For  instance,  a  constitutive  law  to  accurately  describe 
the  behavior  of  snow  under  multi-axial  deformations  needs  to  be  developed.  Such 
a  result  could  find  application  to  a  number  of  problems,  including  vehicle  grade- 
ability,  penetration  mechanics,  shear  wave  propagation,  to  name  a  few.  In 
addition,  there  is  still  a  need  for  more  data,  so  continued  shockwave  testing 
and  high  rate  testing  should  also  be  encouraged  when  possible. 
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